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Foreword 


This  is  part  I of  the  final  technical  report  for  Contract  No.  N62269- 
76-C-0378,  which  is  sponsored  by  the  Naval  Air  Development  Center,  Warminster, 
Fa.  The  work  was  performed  during  the  period  of  July  1,  1976  through 
December  30,  1977.  Hr.  Lee  W.  Gause  was  the  contract  monitor. 

The  contracted  atudy  is  under  the  title  "Certification  of  Composite 
Aircraft  Structures  under  Impact,  Fatigue  and  Environmental  Conditions" ; 
parts  I and  II  of  the  s’-udy  are  under  the  supervision  of  Dr.  P.C.  Chou,  while 
part  III  is  under  Dr.  A.S.D.  Wang,  both  of  Drexel  University. 

This  report  concerns  part  I of  the  contract,  low  apeed  impact  of  plates 
of  composite  materials.  It  is  a self-contained  report,  including  definitions 
of  all  nomenclature  used,  and  its  own  introduction  and  conclusions. 

The  authors  would  like  to  thank  Dr.  Edward  J.  He Quillen,  Dr.  James  l>. 
Huang  and  Hr.  Lee  W.  Gause  for  the  frequent  technical  discussions.  They 
would  also  like  to  thank  Mr.  Frank  Patota  and  Hr.  George  Chou  who  helped 
conducting  part  of  the  experiments. 
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I,  INTRODUCTION 

Aerospace  structures  frequently  undergo  impacts  by  blunt  objects,  in- 
cluding dropped  tools,  hail,  and  runway  stones  impinging  on  exposed  aircraft 
components,  and  foreign  objects  entering  jet  engines.  The  failure  of  a 
component  subject  to  such  an  Impact  is  often  of  the  structural  type,  rather 
than  due  to  local  penetration  or  indentation.  A number  of  numerical  techniques 
(e.g.,  finite-element  method,  lumped-parameter  models)  are  available  for  calcu- 
lating structural  response  to  Impact,  but  these  are  typically  conq> Heated, 
time-consuming,  and  usable  only  on  a problem-by-problem  basis.  Designers  need 
a quick,  convenient,  widely  applicable  method  for  this  purpose. 

In  a previous  report  (1],  & method  was  developed  for  constructing  a de- 
sign curve  which  predicts  the  response  of  a given  type  of  structure  to  impact 
loading.  This  curve  givea  the  maximum  strain  in  the  structure,  which  may  have 
various  dimensions  and  material  properties,  due  to  impacts  Involving  a certain 
range  of  Impact  masses  and  velocities.  An  example  of  the  bending  response  of 
a simply  supported  beam  under  central  impact  was  presented  in  detail.  Doth 
experimental  results  and  numerical  calculations  Involving  several  solutions 
of  beam  impact  were  used  in  establishing  the  design  curve.  The  impact  cases 
studied  in  UJ  were  limited  to  large  impactor  mass,  where  the  lmpactor  has 
mass  roughly  equal  to  or  greater  than  that  of  the  beam. 

In  tba  present  report,  the  design-curve  approach  la  extended  to  anisotropic 
plates,  end  impacts  by  small  lmpactors.  Shear  failures  are  also  studied.  The 
elementary  model  of  Impact  on  beams  presented  in  ( 1]  is  generalized  to  embrace 
all  structures.  This  modal  la  then  applied  to  the  cases  of  clamped  and 
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simply  supported  anisotropic  rectangular  plates  and  to  predicting  impact 
failure  due  to  shear  effects  in  both  beams  and  plates.  Design  curves  are 
also  developed  for  predicting  the  response  of  beams  and  plates  to  impact 
by  small  impactors.  In  each  instance,  both  analytical  methods  and  impact 
experiments  are  employed  in  generating  the  design  curves. 

In  dealing  with  low-velocity  impact  problems,  it  is  convenient  to  divide 
them  into  two  domains,  based  on  the  ratio  M of  the  mass  of  the  structure  to 
the  mass  of  the  impactor  (see  Fig.  1).  The  motion  of  the  impactor  and  the 
response  of  the  structure  for  a large  impactor  differ  greatly  from  those  for  a 
small  impactor.  The  two  domains  will  be  treated  separately;  the  assumptions 
and  final  design  curves  are  also  different  for  the  two  domains. 

In  impacts  where  the  mass  ratio  H is  small  (large  Impactor),  it  has 
been  observed  that  multiple  collisions  occur  and  that  the  general  motion  of 
the  impactor  follows  the  path  of  the  structure  at  the  impact  point  (Fig.  la). 
The  final  rebound  of  the  Impactor  takes  place  after  t.he  structure  has  reached 
its  maximum  deflection.  Further,  the  entire  event  is  generally  more  pro- 
longed than  the  fundamental  period  of  structural  vibration. 

On  the  other  hand,  when  M is  large  (small  impactor),  only  a single, 
sudden  collision  occurs,  after  which  the  impactor  rebounds  and  the  structure 
continues  to  vibrate  freely,  reaching  its  maximum  deflection  at  a later  time 
(Fig.  lb).  The  duration  of  actual  contact  is  characteristically  much  shorter 
than  the  fundamental  period  of  vibration.  The  sudden  rebound  of  the  impactor 
is  due  chiefly  to  the  elastic  resistance  of  the  structure  and  the  Impactor 
to  local  Indentation,  l.e.,  contact  effects  between  the  two  bodies. 

2 


Small  Imp actor 
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In  [1],  six  analytical  models  of  impact  were  examined,  including  two 
classical  one-degree-of-freedom  models;  a two-degrees-of-freedom  model  which 
accounts  for  contact  effects;  a solution  by  Clebsch  [2],  which  assumes  plastic 
impact  and  treats  the  beam  with  attached  impactor  as  a free-vibration  problem; 
a modified  Clebsch  solution  due  to  McQuillen  et  al.  [31;  and  Timoshenko's 
solutiou  [A],  which  couples  the  Impactor  displacement  with  the  beam  deflectiou 
using  Herts's  law  for  contact  deformations.  The  dimensionless  parameters 
which  determine  the  peak  impact  response  according  to  each  solution  method 
were  derived. 

It  was  ahown  in  [1]  that,  for  those  models  which  neglect  contact  effects, 
the  maximum  strain  due  to  impact,  when  properly  normalized  to  e,  la  a function 
of  the  mass  ratio  K only  and  is  Independent  of  any  other  parameters.  These 
solution  methods,  however,  are  generally  not  suited  for  treating  problems 
Involving  small  lmpactora.  Conversely,  the  two  models  that  include  contact 
effects  indicate  a dependence  of  e on  M and  another  parameter,  and  are  useful 
over  the  entire  range  of  M.  For  large  lmpactora,  it  was  demonstrated  that  the 
dependence  on  the  other  parameter  is  weak,  so  that  a single  e vs.  M curve 
approximately  represents  all  impact  situations.  In  addition,  when  data  from 
assorted  impact  experiments  were  plotted  in  coordinates  of  c and  N,  a single 
curve  could  be  drawn  through  all  data  points,  with  about  + 303  variation,  for 
large  impaetors. 
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It  this  report,  several  approaches  are  taken  in  extending  the  design- 
curve  concept  to  impacts  of  other  structures  by  large  impactors  (Section  II) . 
First,  the  one-degree-of-freedom  model  described  in  [1]  is  generalized  to 
encompass  all  structures  and  is  then  applied  to  anisotropic  laminated  rectang- 
ular plates  with  clamped  and  simply  supported  boundary  conditions.  Then,  in 
order  to  identify  additional  significant  parameters,  the  structural  impact 
equations  of  Timoshenko  are  cast  in  a dimensionless  form;  the  influence  of 
each  of  these  parameters  is  appraised  by  a series  of  calculations  using 
Timoshenko's  solution  in  which  the  values  of  the  parameters  are  systematically 
varied.  Ne  , we  discuss  two  series  of  impact  experiments  performed  on  lam- 
inated graphite /epoxy  nlates;  in  one  series,  the  bending  strains  in  the  plates 
were  mes^urud,  and  in  the  ovner,  the  plates  were  repeatedly  Impacted  at  grad- 
ually increasing  velocities  until  failure  occurred.  The  purpose  of  these  ex- 
periments is  twofold:  co  verify  th<*  unalytically  constructed  design  curves,  and 
to  demototvate  now  dee iga  curves  may  aleo  be  constructed  experimentally.  Finally, 
design  curves  for  predicting  failure  due  to  Impact  are  developed,  and 

the  relative  importance  of  shear  and  betting  effects  in  Impact  failure  la  dis- 
cussed. 

For  eaoea  in  which  the  mass  of  the  impactor  is  sjall  with  respect  to  the 
mass  of  the  structure  (largo  H),  the  ' 'en deuce  of  the  generalized  strain  e 
ou  the  mass  ratio  M la  not  exclusive;  <*n„ ther  dimensionless  parameter  whi^h 
Involves  the  geometric  and  material  proparties  governing  the  contact  effects 
between  the  two  bodies  can  also  be  shown  to  be  significant.  This  contact 
parameter,  denoted  X,  end  the  mass  ratio  H together  determine  tue  duration  of 
contact  between  the  structure  and  Impactor  (relative  co  the  fundamental  parlod 
of  vtbraticc),  which  la  typically  quite  short.  Clearly,  the  duration  of  contact 
Is  importune  in  determining  the  response  of  the  structure. 
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The  most  sophisticated  solution  method  applicable  to  problems  involving 
small  impactors,  the  Timoshenko  solution,  involves  a nonlinear  integral 
equation  which  can  only  be  solved  numerically  on  a problem-by-problem  basis, 
and  is  thus  not  convenient  for  developing  a design  curve.  However,  by  making 
the  simplifying  assumption  that  the  structure  does  not  appreciably  deflect 
during  the  short  period  of  actual  contact,  the  problem  uncouples  into  two 
parts  - the  elastic  Hertzian  impact  of  a sphere  on  a semi-infinite  body,  and 
the  vibration  response  of  the  structure  to  a dynamic  load  (i.e.,  the  contact 
force) . This  leads  to  an  approximate  relationship  (for  each  type  of  structure) 
between  the  parameters  e,  H,  and  X,  which  may  be  further  simplified  to  a 
simple  equation  in  the  mass  ratio  M and  a new  generalized  strain  for  small 
impactors  e*  (Section  III) . By  comparison  with  experimental  data  and  with 
calculations  involving  the  complete  Timoshenko  solution,  the  form  of  this  re- 
lationship is  shown  to  be  useful  as  a design  guide. 

In  Section  IV,  procedures  for  using  the  design  curves  presented  in  this 
report  are  described,  and  a few  examples  are  given.  Also  included  are  methods 
for  constructing  new  design  curves  for  structures  not  treated  in  this  report, 
based  on  either  analytical  tools  or  impact  experiments. 
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II.  DESIGN  CURVE  FOR  IMPACTS  BY  LARGE  IMPACTORS 

For  impact  cases  where  the  mass  of  the  structure  is  small  compared  to 
the  mass  of  the  impactor,  the  structure  may  be  assumed  to  behave  as  a simple 
spring,  or  as  a one-degree-of-freedom  mass-spring  system.  Details  of  this 
approach  as  applied  to  simply  supported  beams  were  presented  in  Ref. [1] ; in 
this  section,  this  impact  model  is  briefly  reviewed,  because  it  will  be  useful 
in  the  development  of  parameters  in  plate  impact  problems.  These  parameters 
are  employed  in  constructing  a design  curve  for  predicting  the  response  of 
plates  to  impacts  by  large  masses.  Finally,  by  examining  data  from  a large 
number  of  experiments  and  by  studying  impact  solutions  using  more  sophisticated 
analytical  techniques,  the  validity  of  the  design  curve  is  established  and  the 
range  of  its  usefulness  is  demonstrated. 

A.  Generalized  One-Degree-of-Freedom  Model 

In  this  model,  the  impactor  and  the  structure  are  considered  attached  to- 
gether after  contact  and  move  together  as  a single  mass  m,  and  the  motion  of 
the  structure  is  governed  by  an  equivalent  spring.  The  Initial  velocity  of 
the  combined  mass,  v^,  may  be  determined  by  two  methods:  One  is  based  on 
the  conservation  of  momentum  (Fig.  2) , the  other  on  the  conservation  of  energy 
(Fig.  3).  In  the  first  . case 

m ■ mj  + ea^  (1) 

where  m^  is  the  impactor  mass,  and  em^  is  an  "equivalent"  mass  of  the  structure. 
The  equivalent  mass  may  be  obtained  by  matching  the  total  kinetic  energy  of  the 
structure  with  the  kinetic  energy  of  the  unknown  equivalent  mass  traveling  at  the 
•velocity  of  the  Impact  point  in  tha  structure,  assuming  that  the  deflection  mode . 
shape  la  the  atatlc  deflection  curve.  ,Foc  simply  supported  beams,  , 
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(a)  beiore  impact  (b)  immediately  after  impact 

Fir,,rc  2.  Qnt-da*r«e-cf-fr cadge  modal  of 
impact,  mc—otva  oooaarvad. 


Figure  jj*  Ona-degvcc-of -freedom  model  of 

«*rgy  ’oooMrvcd*  * 

» 
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Vq  - v/ (1  + eM) . 


e ■ 17/35.  By  equating  the  momentum  of  the  impactor  before  impact,  m^  v, 
with  the  momentum  of  the  combined  mass  m,  we  have 

(2) 

where  M ■ m^/m^.  (3) 

Note  that  according  to  Eq.  (2),  the  energy  before  and  after  Impact  is  not 
2 2 

conserved,  i.e. , rav^  < ny/  . 

In  the  energy  conserved  case,  we  assume  e * 0,  so  that 

m - m„ 


V m V 

vQ  V 


(4) 


As  a result,  the  kinetic  energy  in  the  Impactor  is  conserved,  and  will  be 
entirely  converted  to  strain  energy  in  the  structure.  For  small  values  of 
M,  the  difference  between  these  two  approaches  is  small. 

The  spring  constant  Is  the  force  per  unit  deflection,  with  respect 
to  a force  acting  at  the  impact  point  in  the  direction  of  the  impactor 
velocity,  or 

P * K1  W1 

After  acquiring  the  initial  velocity  Vq,  the  mtss-sprlng  system  is  assumed 
to  perform  a free  vibration.  The  maximum  deflection  is  then 

vl,max  “ ^ 

In  modeling  the  structure  by  a mass-spring  system,  it  is  implied  that  only 
the  first  mode  of  vibration  is  retained  and  that  its  mode  shape  is  the  same 
as  the  static  deflection  distribution  under  a concentrated  force. 

Next,  va  shall  consider  the  most  critical  strain  Cj  in  the  structure  which 
occurs  at  a known  point  2.  Assume  that  a proportional  relation  between  tj  end 


can  be  found. 
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where  can  be  determined  from  the  static  deflection  distribution,  or  from 
experimental  static  measurements.  Combining  Eqa.  (5)  and  (6),  we  obtain  the 
maximum  strain  as 


or,  combining  Eqs.  (1),  (2), 


and  (7), 


max 


v if  P2  ~~ 

d12  y<i+eM)Ki 


(7) 


(8) 


1.  Simply  Supported  Beam 

The  next  step  of  defining  a generalized  strain  £ is  best  Illustrated  by 

considering  a specific  type  of  structure.  For  a simply  supported  beam  of 

length  L and  depth  h,  impacted  at  middle  span,  we  have 

d12  - L2/6h  (9) 

and  , 

Kx  - 48RI/L3  . (10) 

Therefore, 


e i 

cmax  % k y4  N(l'feN) 

where  • </e/"d,  the  velocity  of  longitudinal  waves  in  a bar. 
Defining  the  generalized  strain  as 


c 


max  v h 


(U) 


(12) 


we  obtain  finally, 


•w 


H(14«M) 


(13) 


which  is  the  equation  of  the  design  curve  for  simply  supported  beams  impacted 
by  large  masses,  as  given  In  (1J. 


4 


4 


i 
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The  corresponding  expressions  for  d^>  end  e for  a imply  supported 
end  damped  orthotropic  plates  will  be  given  in  the  following  sections. 

2.  Simply  Supported  Plate 

The  impact  design  curve  for  simply  supported,  rectangular, 
orthotropic  plates  may  be  generated  by  applying  the  generalized  one-degree- 
of -freedom  modal.  For  transverse  impacts  at  the  center  of  the  plate,  we  first 
consider  the  corresponding  static  problem  involving  a centrally  applied  con- 
centrated load.  In  this  case,  the  maximum  deflection  occurs  at  the  center  of 
the  plate,  and  the  maximum  strain  on  the  surface  opposite  the  applied  load. 

It  is  shown  in  Appendix  A that  subject  to  a mild  approximation,  the  spring 
constant  of  a plate  so  loaded  may  be  expressed  es 


where 


4ab  fi<n) 


fx(n)  - I l c**" 

4 m-1,3,5  n-1,3,5  mu 


L 2 2 & 2 

C • m + 2m  n n * n V 


the  strain-deflection  cons tent  d^  valeting  the  strain  due  to  bending  in 
the  x-dir action  to  the  maximum  deflection  w^  io 


»*h  f2(«\l 


lab  (fJfjOO 


where 
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f2(n) 


CO 


l 

w»l,3,5 


• 2 

l ?~ 

n*»l,3,5  mn 


(19) 


By  combining  Bqs.  (14)  and  (18)  with  Eq.  (8),  and  assuming  that  e » 0 
(conservation  of  energy) , we  obtain 


*2(n)  hv 


oh  11/2  hv  «2<*> 

'uM  fl(n)  J cpk  [M  f^n)]1^2 


(20) 


If  ve  define  the  generalized  strain  as 


(21) 


where  - *5^/1  is  the  speed  of  flexural  waves  in  the  plate  in  the 
x-direction,  and  k is  the  radius  of  gyration  of  the  plate , then, 

e - (22) 


where 


* r*r 

gjL(n)  • fj^/ZCjCny 


(23) 


This  function  la  plotted  In  Pig..  4,  Similarly,  we  have 

tiCn"1) 


(24) 


Equations  (22)  and  (24)  may  be  plotted  as  impact  design  curves,  aa  shown 

in  Fig.  6.  Mote  that  in  this  esse,  the  generalized  strains  e and  e 

x y 

are  functions  not  only  of  the  structure-to-iapactor  mass  ratio  M but  also 
of  the  aa pec t-orthot ropy  ratio  n. 


3.  Clamped  Plate 

Ve  may  develop  the  design  curve  for  clanped-edge,  rectangular, 
orthotropic  plates  in  a similar  manner,  the  spring  constant  of  such  a plats 
with  r aspect  to  a central  point  load  may  ba  expressed  approximately  (see 


4 


1 


/ 
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Appendix  B)  as 


(25) 

Che  strain  in  Che  x-direction  with 

(26) 

B. 

By  conbinlng  Eqs.  (25)  and  (26)  with  Eq.  (8),  and  assuming  Chat  e ■ 0 
(conservation  of  energy),  and  defining  Che  generalised  strain  according  to  Eq.  (21) 


ve  obtain 

C m 

g2(n) 

(27) 

X 

where 

gjOO  - 

f (n)//f  in) 
n 3 

(28) 

b u 


Kl“ 


11 


a f3(n) 


The  strain-deflection  constant  d^  relating 


the  maximum  deflection  is 


a2  f3(n) 
dl2  " h f, (n) 


where  £3(11)  and  f^(n)  are  given  in  Appendix 


This  function  is  also  plotted  in  Fig,  4.  Similarly,  ve  have 

- «2(n’1> 


y & 

Equations  (27)  and  (29)  represent  the  impact  design  curves  for  a clamped 


(29) 


rectangular  plate*  Again,  the  generalised  strains  e and  c are  functions 

* y 

of  both  the  mass  ratio  M and  the  aspect-orthotropy  ratio  n. 

4.  Comprehensive  Impact  Strain  Curve 

Examination  of  the  design-curve  equations  developed  above,  Eqs.  (22), 
(24),  (27),  and  (29),  shows  that  a different  ve.  H curve  must  be  prepared 
not  only  for  each  boundary  condition  but  also  for  each  value  of  tbs  aspect- 
orthotropy  parameter  n*  These  curves  can  be  combined  into  a single,  com- 
prehensive design  curve  if  ve  Introduce  e new  mans  ion  lass  generalised 
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Figure  Us  function*  of  aspect-crthotrcpy  parameter  for  constructing  iH?>act  design 
carves  for  damped  twi  simply  supported  rectangular  ortbotropic  plates* 
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I 


t 


t 


strain,  e,  defined  as 


_ e 


e - 


. i - 1.2 


(30) 


where  g^(n)  appllos  to  simply  supported  plates  and  is  given  by  Eq.  (23)  and 
gjCh)  applies  to  damped  plates  and  is  given  by  Eq.  (28).  Thus,  the 
equation  for  the  design  curve  in  c vs.  M coordinates  becomes 

T - if 1/2  (31) 

All  impact  data,  regardless  of  the  boundary  condition  (simply  supported  or  clamped) 
or  the  value  of  the  parameter  n may  be  presented  on  a single  curve.  This 
curve  is  shown  in  Fig.  5,  in  which  are  Included  data  from  numerous  impact 
experiments  which  are  discussed  In  detail  in  a later  section. 


5.  Critical  Impact  Velocity  Curve 

The  design  curves  in  terms  of  or  e vs  M give  the  maximum  strain  in  the 
plate  for  a given  impact  situation.  This  strain  may  be  compared  to  the 
ultimate  failure  strain  of  the  plate  (assuming  a naxlmum~strain  failure 
criterion)  to  determine  whether  the  plete  fella  as  a result  of  a given  impact. 
However,  by  Introducing  a new  parameter,  the  dimensionless  Impact  velocity, 
we  can  construct  an  alternate  fora  of  the  design  curve  which  can  be  directly 
used  to  determine  whether  e given  impact  causes  failure  of  the  plete.  Again, 
the  curve  may  be  developed  either  by  analysis  or  by  experiment.  To  do  this. 


e dimensionless  impact  velocity,  v,  is  defined  as 

vh 


v - 


Sk 


(32) 


where  v is  the  impact  velocity,  h is  the  plete  thickness,  is  the  speed  of 
flexurel  waves  In  the  plate,  end  k in  the  radius  of  gyration  of  die  plate. 


IS 


Design  Curve  for  Impact  of  Plate  by  Large  Impactor 


Figure  5.  Comprehensive  impact  strain  curve  for  simply  supported  and 
clamped  rectangular  plates. 
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Then,  by  Eq's  (21),  (22)  and  (32)  the  critical  dimensionless  velocity  v . 

c 

(i.e.,  the  lowest  value  of  v ft  which  failure  occurs)  is  related  to  the  failure 
8 train  by 

\ * £f  <33) 

In  v vs.  M coordinates,  the  curve  v - v£(M)  divides  the  plane  into  a safe 
region  (below  the  curve)  in  which  failures  due  to  impact  do  not  occur  (according 
to  the  assumed  maximum-strain  failure  criterion)  and  an  unsafe  region  (above 
the  curve)  in  which  failures  do  occur.  That  is,  if  the  point  (v,M)  corresponding 
to  a given  impact  situation  falls  above  the  curve,  then  failure  occurs;  If 
the  point  falls  below,  no  failure  occurs.  Thus,  Eq.  (33)  is  useful  as  a de- 
sign curve  equation.  A typical  Vfi  vs.  M curve  is  shown  in  big.  12. 

Further,  if  we  define  a new  dimensionless  Impact  velocity  v such. that 

v - g^n)  v 

then  the  equation  of  the  alternate  design  curve,  in  terms  of  v vs.  M,  la 

v - e.  »H  (34) 

c f 

Mu  wra 

This  single  curve,  v - vc(M),  may  be  used  to  predict  impact-induced  failure 
for  all  place  structures,  regardless  of  the  boundary  conditions  or  the  value 
of  the  aapect-orthotropy  ratio  n (provided  these  structures  have  the  same 
value  of  failure  strain  ef ) . Note  that  the  value  of  efc.  can  be  taken  as  the 
nominal  value  (say,  the  static  ultimate  strain)  of  the  given  material  or  can 

am 

be  determined  by  performing  e tingle  impect-to-feilure  teat.  A typical  v 

c 

vs.  N curve  la  shown  In  Fig.  13. 


B.  Timoshenko  Solution  of  Transverse  Plate  1 


Timoshenko's  approach  for  solving  transverse  impact  problem*  on  simple 
betas  by  coupling  Herts's  lew  of  contact  with  the  Euler  beam  aquation  has  bean 
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extended  to  the  case  of  a simply  supported  rectangular  plate  by  Eringen  [5] 
and  store  recently  to  the  case  of  a simply  supported  anisotropic  laminated 
plate  by  Sun  and  Chattop&dhyay  [6].  In  this  section  we  review  the  salient 
points  of  the  latter  solution  and  show  how  it  may  be  used  to  develop  an  impact 
design  curve  for  anisotropic  plates.  In  addition , a normalized  form  of  the 
equations  corresponding  to  a special  case  of  this  solution  is  derived,  so 
that  the  dimensionless  parameters  governing  the  peak  strain  response  may  be 
identified.  Finally,  a series  of  numerical  computations  is  performed,  in 
which  each  of  the  derived  parameters  is  systematically  varied  in  order  to 
demonstrate  the  significance  of  each  parameter  in  describing  the  response  to 
a given  impact. 

1,  Description  of  Solution  Method 

The  anisotropic  plate  equations  of  Whitney  and  Pagano  [7]  are  used 
to  predict  the  plate  motion.  Neglecting  effects  of  rotatory  inertia,  the 
deflection  of  a symmetric  cross-ply  laminate!  plate  due  to  a centrally  applied 
force  F(t)  is 

grfu-2 

v(x,y,t)  ■~n  f F(x)sinu  (t-i)dt 

■i » n l V h m 

. «in  S2Z  j « 1,3, 5, •••»•*  (35) 

where  m^,  a,  end  b ere  the  mass  end  plener  dimensions  of  the  plate,  and 
uan  *r*  frequencies  dependant  on  a end  tv  end  the  properties  of  the 

plai-u 

£*"*''e  lev  of  oontect  i*  assumed  to  hold 

f • kj  a3'2  (36) 
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where  k2  Is  a constant  and  a is  the  indentation  of  the  impactor  relative 


to  the  plate  surface,  or 


/ft  b \ 
a - w2  - w (j,  t) 


U /) 


for  central  impact , where  w^  is  the  impactor  displacement. 

Newton's  law  applied  to  the  impactor  is 

1 f* 

w-  ■ vt I F dt  dt 

2 m2  ■'0  •'0 

These  four  equations  (35-38)  may  be  combined  into  a single  nonlinear  integral 
equation  in  terms  of  tue  contact  force  F between  the  plate  and  the  impactor. 


.38) 


1 2/3 


,k2.; 


■ vt 


— fC  P F dt  dt 
°2  JO  Jo 


-5“  X _ X < LF(T)  8in  U«n(t“T>dT 
1 m«l,3,5  n-1,3,5  mn  J0 


(39) 


which  is  solved  numerically  by  applying  the  small-increment  method  suggested 
by  Timoshenko 1 in  which  the  contact  force  F is  assumed  to  be  constant  or 
linearly  varying  during  any  time  increment  At.  Expanding  the  above  Integrals 

ah 

to  calculate  the  force  F^  during  the  1 time  interval,  we  obtain 


■i 


2/3 


At2  1 A 1 

- v.At  - E D.  ...  F.  - f-  I B.  ...  F. 
I *2  j*l  i-i+1  J J»1  J 


(AO) 


where 


Bi-j+l  ml  1 “S  {cos(V(1-3)At) 


m n u. 


- coslw^d-J+DAtl},  m,n  • 1,3, 5, 


If  the  contact  force  is  approximated  as  a piecewise  linear  continuous 

th 

function  of  time,  with  an  average  value  of  F^  during  the  1 time  step,  then 
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+ (l-3)(F3-F2+F1)+***+(i-j)<Fj-Fj_1+Fj_2 + Fj)+***+ 


wi-r’i-2+Fi-3'"  - Fi» + 3(ErFi-i+Fi-2'"4Fi) 


For  computing  the  solution  of  Eq.  (40),  Sun  and  Chattopadhyay  have 
suggested  a recursion  method,  but  we  have  found  that  such  a time-saving 

i 

approximation  is  unnecessary.  A listing  of  the  computer  program  for  solving 
this  equation  is  presented  in  Appendix  C. 


2.  Dimensionless  Form  of  Timoshenko  Solutions 
a.  Plate  Impact 

- By  expressing  the  Timoshenko  solution  in  a particular  normalized 
(dimensionless)  form,  the  parameters  governing  the  impact  response  of  plates 
may  be  identified.  Specifically,  it  is  shown  that  the  generalized  strain 
depends  on  three  dimensionless  parameters:  the  mass  ratio  H,  the  aspect- 
orthotropy  ratio  n,  and  also  another  parameter  which  involves  the  geometric 
and  material  properties  governing  the  contact  effects  between  the  plate  and 
the  impactor. 

For  a simply  supported  rectangular  plate  Impacted  at  the  center, 
let  us  consider  the  governing  equation  (39)  and  the  strain  equation, 


NADC  *78259  *60 


2hir 


2 ft 

- — I F(t)  sin  «_(t-T)dT  (41) 

V Jo 

If  we  aasuae  that  the  plate  la  specially  orthotropic  in  beading  (i.e. , 


e « z z 

x ayT  ® n “an  J0 


D26  - 0)  and  further  that  the  flexural  rigidities  are  related  by 


D12  + 2D66  " -jpU  ®22 


(see  Appendix  F)  then  the  natural  frequencies  are 

“an  ’ “U  Am(n) 


(42) 


(43) 


where 


(44) 

Mote  that  this  expression  for  the  natural  frequencies  is  a special  case 
of  that  eaployed  in  the  Sun  and  Chattopadhyay  solution,  in  which  the 
epproxlaatlon  of  Bq.  (42)  is  not  used* 

In  order  to  uoraalise  Bqs.  (39)  end  (41),  the  following  dlaona ionises 
variables  ere  defined, 

**“u‘ 


?.,*•>  .u/Ouvb 

Recalling  the  definition  of  generalised  strain. 


hbSSL 

* hv  "P 


c - c_  k 


(45) 


(21) 


then  Bqs.  (39)  and  (41)  aay  be  written  as 


f f‘  f dt  dt 

/ft  /ft 


w*(Mfir  9 o 


as 


21 


(46) 
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£*  Al+h) 


2 2 f* 

3 £ £ rn-  F(t)sin  A (t-M 

1-hi)3  »n  Am(n)  JQ  an 


where 


and 


M - 


-(¥)’ 


co 

vk 


11 

2 

2 


(47) 


O) 


(48) 


and  t la  a dummy  variable  for  dloenaionless  time  t.  It  can  be  seen  from 
these  normalized  equations  that  the  generalized  strain  depends  only  on 
the  mass  ratio  H,  aspect-orthotropy  ratio  n,  and  the  parameter  Xp;  that  is 


ex  " ex^M,n’ V 

The  parameter  X^  defined  in  Eq.  (48)  may  also  be  expressed  as 


(49) 


(50) 


Thus,  the  parameter  X^  can  be  described  as  the  product  of  several  dimension- 
loss  quantities. 

, / ^ f elastic  forcer  / wave  speed  N 

Xp  - (constant)  x ^ contact  force  J x “impact  velocity/ 

x * length"'^)  * (£uaction  of  «®pect-orthotropy  ratio] 


b.  Beam  Impact 

In  Bef.tll,  the  equations  Cor  the  Timoshenko  solution  of  beam  impact 
were  presented  and  discussed.  These  equations  may  also  be  normalized  by  the 
above  procedure.  If  we  define 

t - w^t 
F - Fl^j/EIv 
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then  the  equations  governing  the  response  of  a beam,  corresponding  to 
Eqs.  (39)  and  (41)  for  the  plate,  may  be  written  as 

CC— 

- 4 5 4 P 5tf)  sin  i2(t-t)dt  (51) 

t 1 iT  J0 
and 


where 


e 


4 I [ F(t)  sin  i2(t-T)dx 
tT  i h 


(52) 


(53) 


which  may  be  impressed  as 


\ - <coMt“t  • 


I deflection  force 
\ contact  force 


■)’ 


/ wave  speed  \ /radius  of  gyration^ 
^impact  velocity/ 


impact  velocity; 


depth  of  beam  J 

It  can  be  seen  from  Eqs.  (51)  and  (52)  that  in  this  case  the  generalised 
•train  c depends  only  on  the  mass  ratio  K and  the  contact  parameter  X^,  or 


e • e(M,l^) 


in  Ref.  (1)  it  was  shown  that  for  large  impact  masses  (M  < -2)  the 
dependence  of  the  generalised  strain  t on  the  contact  parameter  X^  was  only 
slight;  that  is 


e « t(H)  for  small  M. 

Vot  the  case  of  plate  Impact,  the  influence  of  the  contact  parameter  Xp  is 
also  relatively  insignificant  in  determining  the  peak  strain  response  to 
impacts  by  large  masses,  this  will  be  demonstrated  in  the  following  section. 
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3.  Parametric  Study  for  Large  Impactors 

As  can  be  seen  £rom  Eq.  (49),  the  response  of  a specially  orthotropic 
plate  satisfying  the  condition  of  Eq.  (42)  is  described  by  several  param- 
eters: M,  q,  and  Xp.  In  this  section,  we  examine  the  significance  of  each 
of  these  parameters  for  cases  where  the  mass  of  the  impactor  is  large  com- 
pared to  the  mass  of  the  plate.  It  will  be  demonstrated  that,  for  plates 
of  given  aspect-orthotropy  ratio,  n,  the  dependence  of  the  impact  response 
on  parameters  other  than  M is  weak  for  low  values  of  M,  so  that 

“ ex(M,n)  for  small  M. 

To  determine  the  significance  of  the  various  parameters  to  the  impact 
response,  several  plate  Impact  problems  have  been  calculated  using  the 
Timoshenko  solution  method.  Each  problem  of  this  series  is  based  on  an 
impact  situation  with  the  impact  constants  listed  in  Appendix  C.  A para- 
metric study  was  carried  out  by  systematically  varying  the  values  of  the 
impact  velocity  v,  the  Hertzian  contact  stiffness  k^,  the  plate  flexural 
stiffness  matrix  , and  the  impactor  mass  oy  A summary  of  this  study, 
along  with  the  confuted  values  of  the  generalized  strains  ex  and  ey  and  the 
contact  parameter  1^,  is  presented  in  Table  I. 
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TABLE  I 

PARAMETRIC  STUDY 

PLATE  IMPACT  CASES  CALCULATED  BY  THE  TIMOSHENKO  SOLUTION  (n-0.459) 


Mass  Impact 

Ratio,  Velocity, 
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Experiment  - Gr/Ep  plate 
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Figure  6.  Generalised  etrein  curves  for  elaply  supported 

rectangular  orthotropic  plates  (nrO.459)  subjected 
to  central  transverse  impact. 
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Note  that  for  all  of  these  calculations  the  value  of  the  aspect-orthotropy 

ratio  n is  the  sane  (n  ■ 0.459),  An  Inspection  of  this  table  Indicates  that, 

for  any  given  value  of  the  mass  ratio  M,  the  maximum  difference  between  any 

two  values  of  e or  of  e is  only  12%,  whereas  the  largest  value  of  A is  22.6 
x y p 

tines  the  smallest.  It  nay  be  concluded,  therefore,  that  the  dependence  of 
the  generalized  strain  or  on  the  contact  parameter  A^  (or  on  any  other 
parameters)  is  weak,  and  that  for  constant  n a single  ex  vs.  M curve  gives 
an  acceptable  representation  of  all  Impact  cases.  Furthermore,  since  the 
powerful  Timoshenko  solution,  which  includes  the  effects  of  contact  behavior, 
indicates  that  such  behavior  does  not  significantly  contribute  to  the  impact 
response,  then  the  validity  of  the  one-degree-of-freedom  model,  which  neglects 
contact  effects,  is  thereby  substantiated  for  the  range  of  low  M. 

The  Timoshenko  solution  may  be  used  to  construct  an  impact  design  curve 
by  plotting  calculated  points  on  an  eL  vs.  H graph  (Fig.  6).  These  calcu- 
lated  points  alone  can  be  approximated  by  a single  curve  which  may  be  used 
as  a design  curve.  Note  that  such  a curve  would  not  be  very  different  from 
the  curve  corresponding  to  the  one-degree-of-freedom  model  discussed  in  a 
previous  section*  Also  included  in  Fig.  6 are  results  of  a series  of  impact 
experiments  performed  on  the  graphite-epoxy  laminated  plate  described  in 
Appendix  C [ 8 J. 

C.  Impact  Experiments 

Severe!  series  of  impact  experiments  were  performed  on  simply  supported 
end  clamped-edgad  platen,  including  a train-measurement  experiments  and  impact- 
to-f  allure  testa.  The  plate  specimens,  made  of  aluminum  or  graphite/epoxy, 
were  impacted  by  blunt  steel  projectiles  using  a drop- test  apparatus . In  the 
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strain-measurement  experiments,  the  strain  in  the  plate  was  recorded  using 
strain  gages  and  an  oscilloscope.  In  the  impaet-to-failure  tests,  each  plate 
was  impacted  at  gradually  increasing  impact  velocities  until  it  failed.  In 
this  section,  these  experiments  are  described  in  detail,  and  the  results  are 
compared  with  the  theoretically  derived  design  curves. 

1.  Specimens 

The  impact  specimens  included  both  aluminum  and  laminated  graphite/ epoxy 
plates.  Three  6061-T6  aluminum  plates  were  used  with  both  simply  supported 
and  clamped  boundary  conditions.  The  dimensions  and  mass  of  these  plates 
are:  (1)  125x125x6.35  mm,  0.282  kg;  (2)  250x125x6.35  mm,  0.561  kg;  and  (3) 
375x125x6.35  mm,  0.843  kg.  These  dimensions  and  masses  include  only  the  area  be- 
tween the  supports  of  the  plates;  an  additional  10-mm  margin  was  left  along  each 
edge  to  support  the  plates » 

Graphite-epoxy  specimens,  fabricated  from  Hercules,  Inc.,  type  AS/35Q1-6 
pre- impregnated  tape  at  various  lay-ups,  were  used  as  simply  supported  and 
clamped  plates.  The  dimensions,  mass,  and  lay-up  of  each  plate  are  listed  in 
Table  IX.  These  dimensions  and  masses  again  Include  only  the  area  between  the  sup- 
ports of  tha  plates;  a 10- aa  margin  protruded  beyond  the  supports  around  each  plate 

Elastic  properties  for  the  composite  plates  were  celculated  from  lamina- 
tion theory,  using  the.  following  properties  for  each  layer: 


Eu  - 17.7xl06  pal  - 1.220xlOU  H/*: 
S22 


- 1, 3x10*  pel  - 8.96xl09  H/m2 


12 


V23  * °“2 


9 


C12  " °13  “ °* 55x10  J*1  " 3.7'JxlO'  N/» 
G23  - Q.54X106  pal  - 3.72X109  H/m2 


TABLE  II.  CRAPHITE/EPOXY  PLATE  SPECIMENS 
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TABLE  III.  SHEAR  AND  FLEXURAL  STIFFNESSES  OF 
GRAPHITE/ EPOXY  PLATE  SPECIMENS 


Flexural 
Stiffuesses'  D 
(N/m) 


Shear 
Stiffnesaea*  A 
(N/m) 


963. 


-73.3 


-73.3 


748. 


J 

F-series 

H-series! 

19750. 

9970. 

5080. 

4710. 

8100. 

13300. 

-26.7 

-529. 

-26.7 

-267. 

5380. 

4960. 

154000. 


147000. 


155000. 


147000. 
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The  computed  values  of  the  flexural  and  shear  stiffnesses  of  each  plate  are 
summarized  in  Table  HI.  In  all  further  calculations,  the  and  D2q  terms 
are  neglected;  that  is,  the  entire  laminate  is  treated  as  specially  ortho tropic. 
Calculated  values  of  the  parameter  n and  the  functions  g^(n)  and  g^Cn) 
necessary  to  construct  the  design  curves  are  listed  in  Table  XV  for  each 
graphite  epoxy  plate. 

2.  Strain-Measurement  Impact  Experiments 

The  plates  were  centrally  impacted  by  blunt  (6.35  , 12.7,  or  25.4  mm 

contact  radius)  steel  cylinders  of  various  masses  at  several  velocities  using 

f*  drop  test  apparatus.  The  strains,  e and  e were  measured,  in  directions 

m y 

parallel  to  the  plate  edges,  directly  opposite  the  impact  point,  using  a 
Micro-Measurements,  Inc.,  type 84-13-062TT-120 metal-foil  0*-90°  rcsette 
strain  gage,  and  recorded  on  a Tektronix  model  565  oscilloscope.  Impact 
velocities  mere  calculated  from  the  drop  height.  > 

For  impacts  on  clamped  pistes,  the  plates  mere  clamped  along  all 
edges  hy  a frame  of  rectangular  steel  bars  (Fig,  7).  The  entire  assembly  was 
fastened  to  a large  lead  plots  and  rested  on  s concrete  floor. 

In  the  experiments  performed  on  simply  supported  aluminum  plates,  each  plat* 
v*a  rested  on  top  of  a frame  of  rectanguler-erose -section  steel  bars.  For 
the  composite  plates,  however,  the  clasping  device  mentioned  above  was  modi- 
fied to  simulate  s simply  supported  boundary.  Grooves  to  accommodate  3.175  am 
(1/8")  rods  wars  machined  in  each  member  of  the  clamping  device  such  that, 
when  the  clomping  davit*  was  assembled,  the  specimen  plats  would  be  resting 
on  the  rods. 

Experimental  results  of  the  strain-measurement  tests  on  both  clamped 
end  simply  supported  plates  are  presented  in  Appendix  E in  both  tabular  and 
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graphical  (e  vs.  M)  form.  Also  included  in  Appendix  E are  some  typical  strain 
vs.. time  oscilloscope  traces.  Figures  8 and  9 are  typical  examples  of  the  e vs. 

M curves  for  the  clamped  plates  and  the  simply  supported  plates,  respectively. 

In  each  of  these  figures  a dashed  line  has  been  fitted  by  eye  through  the 
experimental  points.  The  maximum  percentage  difference  (variation)  between 
this  line  and  the  experimental  data  is  also  shown  on  each  figure.  For  both 
the  clamped  and  the  simply  supported  cases,  the  largest  value  of  this  variation 
is  36Z. 

Experimental  data  from  specimens  having  approximately  the  same  value  of 
the  parameter  n have  been  presented  on  the  same  design  curve  except  the  three 
smallest  sited  plates  (B1,F1,  and  HI) . Values  of  e for  these  three  plates 
are  considerably  smaller  than  those  for  larger  plates  with  the  same  parameter 
n.  It  is  believed  that,  due  to  the  small  size  of  these  plates,  the  boundary 
constraint  has  more  influence,  and  the  deflection  mode  is  different  from  the  other 
plates.  Perhaps  thiok-plate,  damping,  mod  shear  effects  are  more  pronounced. 

Data  from  these  plates,  although  presented  on  the  i vs.  M plot  of  Figs.  10  and 
11,  were  not  used  in  constructing  the  experimental  design  curve  line  nor  in 
calculating  the  variation. 

Also  shown  In  Figs.  8 and  9,  and  in  each  design  curve  in  Appendix  E,  is 
the  theoretical  (ons-dsgree-of -freedom)  design  curve.  Comparison  of  this 
curve  end  the  one  from  experimental  data  shows  so  average  discrepancy  of  about 
40Z.  Note  that  ths  theoretical  curve  is  always  conservative  (l*e.,  predicts 
e higher  plete  strain)  over  the  range  of  experimantal  results.  However,  the 
elope  of  the  line  drewn  through  the  experimental  data  la  in  all  cases  roughly 
the  same  aa  the  slope  of  ths  theoretical  curve. 
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Figure  11. 


Strain  Curve  for  Simply  Supported  Plates 
Impacted  by  Large  hasses 
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As  prsvioualy  discussed,  ell  inpact  data  can  be  plotted  on  a single 
curve,  7 vs.  M.  This  is  done  for  inpact  data  on  clasped  plates  in  Fig.  10 
and  on  sinply  supported  plates  in  Fig.  11.  Note  that,  since  e is  independent 
of  both  the  parameter  n and  the  boundary  conditions,  both  cases  could  have  been 
plotted  on  a single  figure.  Because  of  the  overlapping  of  experimental  data, 
the  clanped  and  simply  supported  cases,  were  presented  separately  for  clarity. 

3.  Iapact-to-Failure  Experiments 

In  addition  to  the  strain  measurement  impact  experiments,  a series  of 
tests  were  conducted  in  which  laminated  composite  plates  were  impacted  to 
failure.  In  each  of  these  tests,  a plate  was  repeatedly  impacted  by  the 
same  projectile  using  a drop-weight  apparatus  at  gradually  increasing  veloci- 
ties (drop  heights)  until  failure  of  the  specimen  was  detected.  The  specimen 
was  examined  for  failure  after  each  Inpact  both  visually  and  ultrasonically , 
using  a hand  held  pulse- echo  transducer  and  reflectoacope. 

The  avaraga  between  the  highest  Impact  velocity,  Vj,  for  which  a particular 
specimen  did  not  fail,  and  the  lowest  velocity,  vf , for  which  any  failure  was 
datactad,  is  regarded  as  the  critical  or  failure  velocity,  v , or 


V • 

c 


vd  + ,f 


This  velocity  is  used  in  characterising  the  impact  resistance  of  each  specimen. 

The  results  of  these  experiments  are  presented  in  both  graphical  and 

tabular  form  in  Appendix  B.  Figure  12  is  a typical  example  of  the  v vs.  M 

c 

curves  for  the  clamped  plate.  Again,  a dashed  line  has  bean  drawn  to  indicate 
the  design  curve  based  only  on  experimental  results.  A a can  be  seen,  there 
is  good  agreement  between  the  experimental  data  and  the  theoretical  one- 
degree-of-£readom  curve. 

Again,  a design  curve  independent  of  the  perimeter  n and  the  boundary 
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condition*  (v  v*.  M)  may  be  used  to  present  lapse t date.  However , from 
Eq.  (34),  it  is  apparent  that  only  those  plates  with  the  sane  value  of  ef 
should  be  presented  on  this  curve.  For  the  composite  plates,  therefore,  a 
v vs.  M curve  should  be  prepared  for  each  different  lay-up.  Figure  13  is 
such  a plot  for  the  experiments  performed  on  B-eeries  plates. 

After  damage  wan  detected,  several  plates  ware  ultrasonic ally  C-s canned. 
Transducer  output  was  passed  through  an  analog  to  digital  converter.  Data 
was  than  processed  by  S*oitrl<*?  eraaefona  techni^u&s  using  the  first  five  wave- 
forms* Typical  transducer  response  ie  shown  in  Fig.  14  for  an  undamaged 
area  of  a plats.  The  top-eurface  end  bo t ton-surface  echoes  ara  clearly  visible 
with  little  evidence  of  internal  raflactlon.  Fijures  15  and  16  show  trans- 
ducar  raaponsa  over  damagad  areas  of  two  diffarent  spacimana.  Figure  15 
shows  a damagad  r agios  art ending  from  about  one-quaxtar  of  tha  plate  thick- 
ness below  the  surf  see  to  the  midplane  of  the  piste.  This  type  of  damage 
has  baan  obaarved  by  usa  of  the  hend-held  transducer  in  epproalmetely  751  of 
the  failure  teste.  It  la  believed  Chet  this  mode  of  failure  ie  de lamination. 
Figurs  16  la  repraaantatlre  of  tha  raesinlng  251  of  tha  failure  testa.  As 
can  ba  seaa,  thara  is  little  Indication  of  damaga  through  tha  thickness  of 
tha  pints.  This,  along  with  tha  abaanca  of  a bottom-surface  echo  indicates 
that  tha  plate  bee  probably  ben  damagad  near,  tha  top  surface.  It  is  believed 
that  this  mode  of  failure  results  from  fiber  breakage,  due  to  beading  end/or 


contact  effects. 


A structure  subjected  to  Impact  may  also  fail  due  to  the  high  level  of 

•bear  a trees  produced.  In  this  section,  a theory  based  on  the  oee-degrae*of- 

ftesfom  impact  model  is  developed  for  predicting  the  maximum  impact-induced 

transverse  sheer  stress  in  simply  supported  beams  end  plates. 
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plate  specimen. 
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1.  Simply  Supported  Beam 

By  extending  the  one-degree-of -freedom,  energy-conserved  model  previously 
discussed,  4 design  curve  for  predicting  the  peak  shear  stress  in  a simply 
supported  beam  subjected  to  transverse  impact  may  be  generated. 

First,  the  maximum  shear  stress  in  the  beam  is  related  to  the  maximum 
deflection  by  the  static  besm  equation.  The  deflection  of  a simply  supported 
beam  subjected  to  a static  load  at  midspan  is 


v(x)  - wx  f (3  - 0<x<L/2 


(54) 


where  is  the  deflection  at  midapan.  The  shear  force  generated  is 


V - BI 


d3v 

to3 


1 


(55) 


For  a beam  of  rectangular  cross-section,  the  maximum  shear  stress  occurs  at 
the  neutral  axis  and  is  given  by 


T - lx 

max  2A 


(56) 


The  maximum  deflection  v 


\mx 


36BI 
AL3  fvx 

may  be  related  to  the  impact  velocity  by 


the  cooservatioo-of-energy  condition,  which  laplies  that  the  initial  kinetic 


energy  of  the  iapactoc  is  entirely  converted  into  banding-etrein  energy  in 
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By  defining  the  dimensionless  (generalized)  shear  stress  as 


t L a A 
max 

vEl 


(58) 


and  combining  Eqs.  (56)  and  (57) , we  obtain  the  equation  for  the  shear-stress 
design  curve. 

7 - 3,/f  (59) 

This  equation  is  plotted  in  t vs.  M coordinates  in  Fig.  17. 

a.  Relative  Importance  of  Shear  and  Bending 

Depending  on  geometry,  strength,  and  elastic  properties,  a particular 
beam  subjected  to  low-velocity  impact  may  experience  failure  Initiated  by 
bending  or  shear  effects.  The  theories  presented  here  for  predicting  bending 
strain  and  shear  stress  due  to  impact  are  both  based  on  the  static  beam 
deflection  curve.  Therefore,  the  relative  Importance  of  bending  and  shear 
effects  msy  be  estimated  by  considering  the  static  relations  between  shear 
and  bending. 

If  a maximum-stress  failure  criterion  is  adopted  for  both  shear  and 
bending,  then  failure  due  to  shear  will  be  more  likely  if,  for  any  impact 
velocity  v,  wa  have 


~ > -2. 


(60) 


‘f  vf 

where  the  subscript  f denotes  values  at  failure.  By  combining  this  relation 
with  Eqs,  (12),  (13,  (5S),  and  (50),  we  obtain 


L Tf 


(61) 


Based  on  the  approximate  theory  presented  hsre,  this  inequality  must  be 
satisfied  for  impact  failure  to  initiate  by  shear  effects;  if  it  is  not 
satisfied,  failure  will  initiate  by  bending. 
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b.  Critical  Impact  Velocity  Curve 

The  relative  importance  of  shear  and  bending  effects  in  initiating 
failure  of  a particular  beam  nay  be  illustrated  by  presenting  the  design 
curve  in  a form  similar  to  that  presented  in  Section  II  for  plate  Impacts. 

In  this  alternate  form,  theoretical  design  curves  and  experimental  data  are 
plotted  in  coordinates  of  a dimensionless  impact  velocity  v or  v and  the  mass 
ratio  M.  Then,  each  critical  velocity  curve  divides  the  plane  into  two  regions: 
above  the  critical  bending  curve  is  the  region  bending  failure  occurs,  and 
below  is  the  no-bending-failure  region.  The  same  is  true  for  the  critical 
shear  curve. 

Specifically,  if  the  dimensionless  impact  velocity  for  the  beam  is  defined 


as 


vh 

V._ 


(62) 


then  the  one-degree-of-freedom  energy-conserved  impact  model,  predicts  that 
failure  will  Initiate  in  the  beam  due  to  bending  for 


2c. 


(63) 


where  ef  Is  the  tensile  failure  strain,  assuming  a maximum-strain  failure 
criterion. 

Similarly,  according  to  the  theory  discussed  above  for  estimating  the 
shear  stress  due  to  Impact,  Eq.  (59),  failure  due  to  shear  occurs  for 


t * 4L 


(64) 


where  is  the  shear  stress  at  failure,  assuming  a maximum-stress  failure 
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Equation*  (63)  and  (64),  plotted  in  v vs.  M eoordinatss,  rapressnt  the  two 
critical  velocity  design  curves.  The  relative  iaportance  of  shear  and  bending 
regarding  failure  of  a particular  structure  is  immediately  apparent  from  such 
a plot  — the  lover  curve  predicts  failure  at  a smaller  Impact  velocity,  indi- 
cating that  the  failure  mode  (shear  or  bending)  associated  with  the  lover 
curve  is  the  critical  one. 

Consider,  for  example,  the  Impact  failure  tests  on  22  graphite-epoxy  beams. 
These  beams  were  fabricated  from  Hercules  AS3501  at  a lay-up  of  [+  45/02/+  4512* 
411  beams  have  a span  to  thickness  ratio  (L/h)  of  48.53.  Data  from  all  the 
impact  tests  (with  or  without  failure)  performed  on  these  beam  specimens  are 
tabulated  in  Appendix  E and  plotted  in  v vs.  M coordinates  in  Fig.  18.  Also 
Included  in  this  figure  are  curves  representing  Eqs.  (63)  and  (64).  The  value 
of  the  bending  failure  strain  used  is  0.0168;  the  value  of  the  shear  failure 
stress  used  is  the  value  of  the  Interlaminar  shear  strength  given  by  the 
material  manufacturer  (Hercules  Product  Data  Sheet  No.  832)  as  18,900  pal 
(130.  MN/m2). 

Note  that  the  Iwer  curve  in  Pig.  18,  which  represents  Eq.  (63),  approxi- 
mately divides  the  experimental  points  according  to  whether  die  points  correspond 
to  failure  or  no- failure  tests.  Also,  the  curve  corresponding  to  shear  failure 
falls  much  higher  on  this  plot  than  the  curve  for  handing  failure;  this  indicates 
that  for  thasa  beam*  the  bending  is  the  mode  in  which  failure  inltietee. 

2.  Simply  Supported  Plate 

A design  curve  for  eat lasting  the  peak  transverse  shear  stress  at  the 
edge  of  e elnply  supported,  rectangular,  orthotropic  plate  subjected  to  central 
lateral  Impact  may  also  be  developed. 

It  must  be  recognised  that  the  edge  shear  la  only  a rough  estimate  of 
the  maximum  shear  stress  occurring  in  the  plate;  such  higher  stresses  may  be 

generated  is  the  immediate  vicinity  of  the  impact  point.  However,  analysis 
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of  the  problem  by  plate  theory  aloneXas  is  done  here)  leads  to  a singularity 
in  the  shear  at  this  point,  due  to  treatment  of  the  impact  load  as  a concen- 
trated force.  A more  sophisticated  analysis  would  focus  'on  the  three-dimen- 
sional stress  field  near  the  impact  location.  Still,  for  design  purposes, 
knowledge  of  the  edge  shear  can  sometimes  be  useful.  If,  for  a particular 
impact  situation,  it  exceeds  the  allowable  shear  stress,  failure  will  certainly 
occur;  on  the  other  hand,  a low  edge  shear  stress  does  not,  of  course, 
guarantee  survival  of  the  structure. 

As  in  the  previous  case  of  the  beam,  the  maximum  edge  shear  stress  in 

the  plate  is  first  related  to  the  maximum  central  deflection  by  asstiming  that 

the  static  relations  hold.  In  terms  of  generalized  plate  forces,  the  maximum 

values  of  the  transverse  shear  stress  components  are 

3Q 


xz  2h 

t .2l 

yz  2h 


(65) 


for  a homogeneous  plate.  For  a specially  orthotropic  plate,  the  generalized 
forces  are  related  to  the  deflection  distribution  by 


-4-4 


32w 


32w  ! 


where 


o.—e-utrv. 

Q L H iJi  + T> 

**  * ax2  22  ay2 


H ■ °12  * 2D66 


(66) 


It  is  again  assumed  that  the  dynamic  deflection  due  to  impact  is  the 
tame  as  the  atatic  deflection  due  to  a central  point  load.  In  this  case,  the 
maximum  edge  thaar  force  occurs  at  the  middle  of  the  edge.  If  it  is  also  assumed 


XX 
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a .3  *i« 


3.3d  yrt  ^ 

2h  b3  f^if1)  3! 


(67) 


(68) 


where 


f5(ti) 


i i ^4- 

art,3,5  n-1,3,5  m*  + n\ 


(69) 


and  where  f^(n)  and  n are  defined  in  Eqe.  (13)  and  (16) , respectively* 

If  we  define  the  dtesnsioalass  sheer  stress  parameters 

2 


Txs"  T, 


V* 
D11Y 
‘ bh 


ya  y*^  DjjjV 

Nx 

then  Bqs.  (67)  and  (68)  nay  be  combined 


(70) 


h the  cooaervatioo-of-enerty 

conditlon«  Eq.  ( 4 ) * to  yield  the  design  curve  equations. 


83(h) 


where 


X**m-xr 

«3(,rl> 
t3«5(n) 


y* 


83(n)  - 


(73) 


Equation  (71)  is  plotted  in  Fi*.  19  for  a few  value*  of  n»  end  the  function 
$3(n)  is  presented  in  Fi*.  20* 
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supported  rectangular  orthotropic  plates 


III.  IMPACTS  BY  SHALL  IMPACTORS 

In  this  section,  approximate  equations  for  predicting  the  peak  strain 
response  of  simply  supported  beams  and  plates  to  impacts  by  small  impactors 
are  derived.  The  accuracy  of  these  equations  is  demonstrated  by  comparison 
with  experimental  results  and  with  the  Timoshenko  solution. 

A.  Simply  Supported  Beam 
1.  Approximate  Analysis 

It  has  been  observed  in  experiments  and  in  calculations  of  the 
Timoshenko  solution  that  for  the  case  of  small  impactors  (mass  retio  K 
greater  than  about  2),  the  Interaction  between  beam  and  impactor  is  a single, 
sudden  blow  of  short  duration  compared  to,  say,  the  fundamental  period  of 
vibration  of  the  beam.  This  Is  quite  different  from  the  large  impactor  case 
where  multiple  Impacts  occur.  By  assuming  that  the  beam  does  not  appreciably 
deflect  during  the  short  period  of  contact  and  by  making  a few  other  approxi- 
mations, ve  may  derive  a relation  vhlcb  la  useful  in  generating  design  curves. 

The  assumption  that  the  beam  does  not  deflect  during  contact  is  equiva- 
lent to  assuming  that  the  impactor  is  much  less  massive  than  the  beam,  which 
may  thus  be  treated  as  & semi-infinite  body.  An  approximate  solution  for  the 
contact  force  during  the  elastic  l&pact  of  a spherical  Impactor  against  a 
fist  semi-infinite  body  has  been  calculated  by  Hunter  [9]  as 

V * [?.  8in  £t  0 < * ! *{l  (74) 

(75) 

(76) 


u 


C » 


where 


5>  - 1.068  v 2 6 

* . o.sna{~^j 
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The  response  of  a simply  supported  beam  to  a centrally  applied  trans- 
verse dynamic  load  of  the  general  form  of  Eq.  (74)  is  derived  in  Appendix  Jl. 
The  maximum  deflection  end  bending  strain  occur  at  midspan  at  a time  about 
one-quarter  of  the  fundamental  period  after  impact.  The  peak  strain  may  be 


expressed  as 


where 


e - “ f(e) 

123  X EX 


*/nv  1 r 23  i tt 

f(6)  - -y  l 2 ~ cos  ~2& 

v i*J.,3,5  & - i a 


1 fundamental  period  of  beam 


contact  force  duration 


Substitution  of  the  approximate  value  of  £ given  in  Eq.  (V6)  into 


Eq.  (79>  yields  tn  approximate  value  for  B. 

,L  2\rhzi/5 

S « (0.09897M  h— ^ J 


This  quantity  may  be  expressed  in  terms  of  the  dimensionless  parameters 
previously  derived  as 


2 v 

6 - (0,1564) 4 r~  ) 
*b  ' 


where  H and  V,  are  defined  in  Eqs.  (3)  and  (53),  respectively. 
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By  substituting  for  P and  B into  Eq.  (77)  according  to  Eqs.  (75)  and 
(76),  and  recalling  the  definition  of  generalized  strain,  Eq.  (12) 
we  obtain 

e - 1.068  it2  f($  ) £ (82) 

This  result  nay  itself  be  used  to  generate  a series  of  design  curves  by 
plotting  e vs.  M for  different  values  of  8 , however,  this  equation  may  be 
further  simplified  as  demonstrated  below. 

Note  that  Eq.  (82)  shows  a strong  dependence  of  the  generalized  strain 
e on  both  the  mass  ratio  M and  the  parameter  8 . Also,  observe  that  for  a 
constant  value  of  8 , e is  inversely  proportional  to  M,  rather  than  vfi 
as  for  the  case  of  large  impactors.  This  suggests  that  for  small  impactora 
the  peak  strain  is  related  to  the  momentum  of  the  impact or,  rather  than  to 
its  kinetic  energy  as  for  large  impactors.  That  is,  large-M  cases  are  governed 
by  the  impulse  exerted  by  the  impactor,  whereas  saall-M  cases  are  character- 
ized by  the  work  done  on  the  beam. 

The  influence  of  the  value  of  the  parameter  8 on  the  generalised 
strain  e can  be  estimated  by  taking  the  limit  of  Eq.  (82)  as  8 becomes  large. 

In  Table  V,  calculated  values  of  f (8  )are  presented  (see  Appendix  D). 
Examination  of  this  table  reveals  that,  as  8 increases,  the  quantity  /f5  f (8) 
approaches  the  value  0.12412.  Observe  also  that,  for  8 > 7,  /(i  f(8)  never 
differs  from  this  value  by  more  than  2X.  Since  roost  practical  problems  satisfy 
the  condition  8 > 7,  we  can  assume  with  little  inaccuracy  that 
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TABLE  V 


CQMPU1  CR-CA1.CULATF.D  VALUES  OK  f(g) 


30.0 

50.0 
oo.o 

100.0 

150.0 

300.0 

300.0 

500.0 


0.06466 

0.0 4 590 

0.03P03 

0.03150 

0,02 ‘/SO 

0.02264 

0.03/55 

0.01338 

0.03  2413 

0.010132 

0.008/7? 

0.007106 

0.003351 


G m) 

0.11200 

0.12145 

0.12342 

0.12356 

0.12432 

0.12402 

0.1240& 

0.12414 

0.12413 

0.12412 

0.1241 3 

0.12412 

0.1241? 
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Substitution  into  Eq.  (82)  yields 


e 1.308 

/r“  M 


(83) 


or,  defining  a new  generalized  strain  for  the  small  impactor  case  as 

e*  - Tjr  <84> 

we  then  have 


e* 


1.308 

M 


(85) 


This  simplified  equation  may  be  used  as  a design  curve  for  the  case  of  small 
impactors  by  plotting  e*  vs.  M.  In  Fig.  21,  such  a curve  is  compared  with 
the  experimental  data  recorded  in  Appendix  E.  In  general  the  experimentally 
determined  values  of  e*  are  roughly  60%  lower  than  those  predicted  by  Eq.  (85). 
Note,  however,  that  the  band  of  experimental  points  on  this  graph  form  a 
relatively  narrow  band;  for  any  value  of  M,  the  veriatlon  in  e*  is  at  most 
+ 25%  (at  about  M ■ 6)  and  much  less  for  the  higher  range  of  M (greater  than 
about  10).  This  variation  in  e*  is  much  smaller  than  the  variation  in  e 
alone,  for  the  same  data.  Also,  the  slope  of  this  band  is  very  close  to 
that  of  the  curve  corresponding  to  Eq.  (85).  This  suggests  that  the  locus 
of  the  experimental  data  plotted  in  these  particular  coordinates  may  Itself 
be  useful  as  a design  curve.  We  may  conclude  that  Eq.  (85),  while  not  accurate, 
is  significant  in  that  it  serves  to  Identify  the  form  of  the  relationship  among 
the  several  pereaetere  governing  beam  response  to  impacts  by  small  masses. 

2*  Timoshenko  Solution 

Also  Included  in  Fig.  21  are  several  points  corresponding  to  calcula- 
tions using  the  Timoshenko  solution;  data  from  these  calculations  are  summar- 
ized in  Table  VI.  Note  that  these  points  form  e very  narrow  band,  the 
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TABLE  VI 


CALCULATIONS  USING  TIMOSHENKO 
SOLUTION  FOR  LARGE  MASS  RATIO  M 


Beam  Mass, 
Dimensions, 

Mass 

Ratio, 

Impact 

Velocity, 

and  Material 

v (m/sec) 

0.389  kg 

27.7 

8.72 

197  x 16  x 16  can 

27.0 

6,37 

tt 

l 13.17 

steel 

13.84 

8.72 

II 

12.19 

•1 

8.72  ■ 

II 

17.43 

I 

II 

8.72  1 

6.86 

5.49 

It 

8.53 

3.54 

— 

II 

9.20 

0.389  kg,  steel 

394  x 6 x 16  nm 

13.84 

8.72 

0.389  kg,  steel 

27.7 

8.72 

394  x 16  x 16  nm 

6.92 

8.72 

0.389  Up,,  steel 

197  x 8 x 32  r.n 

13.84 

8.72 

0.1943  kg,  $tccl 

197  x 8 x 16  rm 

6.92 

8.72 

2.466x10 


2.466x10 


l . 744x10 


2 .41 6x1 Ol 


2.466x10 


Generalized 

Strain, 

e ■ ca2/hv 


0.1807 


2.415xl01U  0.1782  ' 


0.1894 


0.2815 


0.2909 


0.2557 


0.2999 


0.2696 


0.4124 


0.4240 


0.5727 


0.5665 


0.2502 


0.0405 


0.0416 


0.0412 


0.0734 


0.0728 


0,0757 


0.0720 


0.0744 


0.1237 


0.1266 


0.3  on 


0.1984 


0.0322 


0.2703  0.0361 


0.3171  0.3559 


0.2163  0.0789 


0.4824  0.1192 
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Figure  21.  Generalised  strain  c*  va.  Mas  ratio  M for  lapacts 
of  elaply  supported  beaas  by  asall  lapactor, 
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variation  being  less  than  + 5Z;  again,  the  variation  in  e*  for  the  sane 
data  is  ouch  larger.  In  addition,  the  band  of  points  is  ouch  closer  to  the 
curve  corresponding  to  £q.  (85)  than  are  the  experimental  points;  these  points 
lie  only  15  to  45%  below  the  curve  and  appear  to  approach  the  curve  as  H 
increases.  The  closeness  of  these  results  implies  that  the  lack  of  agree- 
ment between  Eq.  (85)  and  the  experimental  data  is  not  only  due  to  the  approx- 
imations used  in  deriving  the  equation.  Apparently,  there  are  additional 
effects  which  significantly  influence  beam  response  to  impacts  by  small  masses 
which  are  not  accounted  for  even  in  the  powerful  Timoshenko  solution.  Several 
researchers  have  modified  the  basic  solution  to  include  some  of  these  effects, 
but  these  more  involved  models  are  too  cumbersome  for  the  development  of  an 
impact  design  curve,  since  for  each  effect  considered,  another  parameter  (or 
possibly  several)  must  be  introduced.  Such  complications  are  surely  justified 
when  one  is  analysing  a specific  i^>act  problem,  but  are  of  little  assistance 
when  one  is  designing  a structure  to  resist  potential  Impacts  due  to  a range 
of  masses,  velocities,  materials,  etc. 

B,  Simply  Supported  Plate 

In  this  section  is  developed  s design  curve  for  predicting  the  response 
of  a simply  supported  piste  to  impact  by  a small  mesa  using  the  same  approach 
as  applied  above  to  a beam, 

1.  Approximate  Analysis 

As  in  the  case  of  the  beam,  we  assume  that  the  plate  does  not  appreciably 
deflect  during  the  short  period  when  the  impactor  is  actually  in  contact  with 
the  plate.  Again,  this  assumption  is  equivalent  to  supposing  that  the  impactor 
Is  much  less  massive  than  the  plate.  Thus,  the  contact  force  is  roughly  given 
by  the  epproxlmtte  solution  for  e sphere  elastically  striking  s flat  Hal- 

iafintts  body,  Kqa.  (74),  (75),  and  (76). 
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The  dynamic  response  of  a simply  supported  plate  to  a central  transverse 
load  F(t)  is  given  by  [ 6 ] 

v.±ll  C-a("^'.2)-/2  6lp aa ,to sa 

m*  " • u>  & d 

1 in  n mn 


II 


F(t)  sin  u (t-T)dx 


(86) 


For  a thin  specially  orthotropic  plate,  the  natural  frequencies  are  given 
by  110] 


(i) 


2 it' 


an  ph 


DU  7 + 2CD12+2D66>  77  + D22  b4 


(87) 


If  we  again  assume  that 


then  we  have 


in  which 


D12  + 2D66  " fu®22 


2 2 
a •»  nn 

“an  “ll  1+q 


2 ,D .,\/2 

“ll  * "2(7^)  (1+f0 


(43) 


(88) 


Substitution  of  the  assumed  fora  of  the  contact  force,  Eq.  (74),  into 
Eq.  (66)  and  integration  yields 


atfn-2 


2 •*»* 

"l  11  a 0 a4  + qn  * 


*4**  SJ£ 


sin 


u ...  * 


• in  w^U-fr/O  cos 


e(i-«  ‘/c  ) 


(89) 


In  the  corresponding  formula  for  the  beam,  all  terms  in  the  aeries  reach 
their  maximal  values  at  erectly  the  seas  time,  due  to  a simple  relationship 
which  exists  among  the  natural  frequencies  (sea  Appendix  D,  Eq.  (01)). 
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In  the  present  case  of  the  plate,  the  terns  In  Eq.  (89)  do  not,  in  general, 
reach  their  peak  values  simultaneously;  however,  we  can  estimate  an  upper 
bound  for  the  deflection  (and  also  the  bending  strain)  by  letting 


sin  to  (t  + tf/O  - 1 

Tnn 


(90) 


for  all  m and  n. 

The  maximum  strain  in  the  x-direction  due  to  bending  occurs  at  the 
center  of  the  plate. 


h 32w 

2 ax2 


x.y,! 
a b 2 


(91) 


so  that  the  upper  bound  on  the  strain  is  given  by 


4_Ph«  1 v r m2  5e. 


II  5- 


irA 


x,Bax  rV.b  (1+n)2  an  Amn  62  - A 2 
11  p nn 


cos 


mn 


28. 


where 


and 


% - *'«u 


mn 


2,2 

a + nn 

1+n 


(92) 


(93) 


(94) 


Substitution  of  the  approximate  value  of  £ given  in  Eq,  (76)  into 

Eq.  (93)  yields  an  approximate  value  for  8 » given  by 

p 1/5 


8 - 0.09897 
i> 


(95) 


This  quantity  may  be  expressed  in  tors*  of  dimensionless  parameters  previously 
derived  as 

1/5 


8 - 0.1564 


H 


\\  (1+n) 


(96) 
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Now,  if  we  substitute  for  P in  Eq.  (92)  according  to  Eq.  (75) , we 
obtain 


x,max 


(97) 


Finally,  if  we  define  the  generalized  strain  for  impact  of  a simply  supported 
plate  by  a small  mass  as 


'x.max  r 11 


T^TTTph 


*x  4(1.068)hv  8p  g(8p,n) 


(98) 


where 


8<Vn)  ’^{{f 


8 


v A 


„ 2 .2 

mn  8 ~ A 

p mn 


cos 


mn 


28 


(99) 


then  the  equation  for  the  design  curve  according  to  this  approximate  solution 


la 

e*  - 1/M  (100) 

Note  that  the  function  g(8p»n),  plotted  In  Fig.  22,  is  practically 

independent  of  the  value  of  8 , at  least  for  realistic  values  of  0 (6_  > 7). 

P P P 

In  Fig.  22a,  g(Bp,n)  is  plotted  over  n for  two  widely  different  values  of 

8 ,7  and  500;  the  two  curves  are  almost  the  same.  In  Fig.  22b,  the  maximum 
P 

and  minimum  values  of  g(8  ,n)  for  8 in  the  range  7 £ 8 < SOG  are  plotted 

P P P 

for  each  value  of  n;  again,  the  two  curves  are  quite  close,  indicating  that 

g(8_,n)  doaa  not  vary  much  with  8 . 

P p 
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Figure  22.  Plot*  of  shoving  weak  dependence  on  0 for  7 < 0 < 500 
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2 . Comparison  with  Experimental  Results 

In  Fig.  23,  the  design  curve  £or  impact  of  plates  by  small  masses  is 
presented  along  with  the  results  from  experiments  performed  at  Drexel.  Two 
graphite-epoxy  plates  which  were  also  used  in  the  large  Impactor  experiments, 

FI  and  F2,  were  impacted  by  aluminum  projectiles.  These  experiments  were 
performed  using  an  air  gun  with  the  projectile  velocity  recorded  by  a photo- 
diode system.  Data  from  these  experiments  are  tabulated  in  Appendix  E.  As 
can  be  seen  in  Fig.  23,  the  agreement  between  the  theoretical  design  curve  and 
the  design  curve  based  on  experimental  results  is  of  the  same  order  as  in  the 
large  impactor  case;  that  is,  the  experimental  curve  is  about  40%  below  the 
theoretical  curve  and  has  approximately  the  same  slope. 

Also  tabulated  in  Appendix  C are  experimental  data  found  in  Schwieger 
(11] . These  experiments  were  performed  on  a square  Duraluminum  (AlCuMgl)  * 
plate,  330  x 530  x 4.97  mm,  centrally  impacted  by  steel  spheres  of  20-  and 
30-  tm  radius.  These  data,  whan  plotted  i;i  e*  vs.  H coordinates,  do  not  compare 
well  With  either  the  Drexel  experiments  or  the' theory.  One  explanation  of  this 
discrepancy  lies  in  the  siae  of  the  plate  tested  by  Schvicger.  In  the  Drexel 
experiments  i the  ratio  of  the  plate  span  to  thickness  (l./h)  has  a value  of  18 
whereas  in  the  experiments  performed  by  Schvioget,  this  ratio  has  a value  of  110. 
Thus,  according  to  Zener  [12],  the  response  of  the  plate  used  by  Schwinger  is 
governed  by  membrane  and  wave  propagation  effects.  Since  these  effects  are 
not  treated  by  the  current  theory,  the  data  from  Schwieger  are  not  included 
on  the  design  curve. 
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IV.  DESIGN  PROCEBUivE 


In  this  section,  procedures  are  recommended  for  using  the  design  curves 
presented  in  this  report  and  for  constructing  new  design  curves  for  structures 
not  treated  here  based  either  on  analytical  calculations  or  on  experimental 
data.  It  la  emphasized  that  these  design  curves  are  limited  to  predicting 
only  the  maximum  structural  response  to  low-velocity  impact.  ' 

A.  Use  of  Design  Curves  in  This  Report 

t 

The  design  curves  in  this  report  predict  the  peak  bending  strain  in  simply 
supported  and  clamped  plates  impacted  by  large  impactors,  and  in  simply  supported 
beams  and  plates  impacted  by  small  impactors.  Also,  in  [1],  simply  supported 
beams  impacted  by  large  masses  were  treated  in  detail.  Methods  for  the  use  of 
each  of  these  design  curves  will  be  discussed  separately  here. 

1.  Simply  Supported  Seam 

Design  curves  for  simply  supported  beams  impacted  by  large  impactors  and 
by  small  impactors  are  repeated  here  for  convenience  (Figs.  24  and  2S,  respec- 
tively). In  each  case,  ve  have  selected  a design  curve  fit  by  eye  to  the 
experimental  results  ae  the  one  ve  will  use. 

The  first  step  in  using  these  curves  Is  to  compute  the  mass  ratio 


aX 

H • rr* 

*2 


icspactor  mass 


(101) 


If  K is  lass  than  2,  the  impact  case  falls  into  the  domain  of  large 
Impactors;  if  H > 2*  sasll  itepaerora. 
e.  Large  Isgactars 

From  Fig.  24,  read  the  value  of  the  generalized  strain  c corre- 
sponding to  the  value  of  M.  then,  using  this  value  c,  compute  the  maximum 
Impact  strain. 
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Figure  25. 


MASS  RATIO, 


Selected  deelgn  curve  for  lapects  of  slaply  supported 
betas  by  Basil  lapse tors. 
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- hv 
Emax  kc 

where  h * beam  depth 

k - radius  of  gyration  of  beam 
cross-section 

v ■ impact  velocity 

c * »/e/p  ■ speed  of  flexural  waves 
in  the  beam 

The  use  of  this  curve  will  be  illustrated  by  an  example. 


(102) 


Example  #1 

3 

For  a simply  supported  steel  (density,  7.9  g/cm  ) beam  of  rectangular 
cross-section  with  dimensions  Lxhitb*  500x25x20  nm  impacted  at  10  m/s 
by  a mass  of  4.0  kg,  we  have 


M 


°2 


L bh 
“2 


(7900)  (.50)  (0.025X0.020) 
4.0 


- 0.494 

Prom  Fig.  24,  the  value  of  e corresponding  to  M ■ 0.494  is  £ ■ 1.20.  Then, 
by  Eq.  (102),  we  have 


e 


max 


(1  20)  <M25lOfiL_ 

/0. 025^  210xl0\1/2 

\7n~  7900*7 


* 0.0081 


b.  Small  Impactora 

From  Fig.  25,  read  the  value  of  the  generalised  strain  e*  corres- 
ponding to  the  value  of  M.  Then  compute  the  maximum  of  impact  strain. 


e^_  ■ c 

max 


ck 


(103) 


where 


6 - (0.09897)  ~ 


and  k2  ■ Hertzian  contact  stiffness. 
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Example  #2 

Consider  the  beam  of  Example  #1  impacted  by  a steel  sphere  of  radius 
20  mm  at  20  m/s.  Then,  by  Eq.  (101),  we  have 

M - i - (7900) (0.50) (0.025) (0.020) 

®2  | * (0.020) 3 (7900) 

■ 7.46 

From  Fig.  25,  the  value  of  e*  corresponding  to  M ■ 7.46  is  e*  • 0,060.  The 
Hertzian  contact  stiffness  is 


so  that 


3ir  Lit  210x10S  J 


2.176xl010  N/m3/2 


0 - (0.09897) 


i 


r - VT/5 


(0.50)*:  20(2.176x10^) 

210x10^  0.025  L 0.2652  J 


7900  Vl2 


- 28.1 


Therefore,  the  maximum  bending  strain  Is 


e - (0.060)  vTs.l 

XmX 


(0.025)  (20) 

if 210xl09  (0,025) 
V 7900 


i 


- .00427 
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2.  Simply  Supported  and  Clamped  Plates 

The  procedures  for  using  the  design  curves  for  rectangular  plates  is 
similar  to  those  outlined  above  for  the  beam.  For  convenience,  we  have  re- 
peated here  the  design  curves  for  simply  supported  and  clamped  plates  Impacted 
by  large  impactors  (Figs.  26  and  27),  and  for  simply  supported  plates  impacted 
by  small  Impactors  (Fig.  28).  The  first  step  in  using  any  of  these  curves  is 
to  compute  the  mass  ratio  M according  to  Eq.  (101)  and  the  aspect-orthotropy 
ratio, 


where  a and  b are  the  plate  dimensions  in  the  x and  y directions.  For  the 
plates  treated  here,  impact  cases  in  which  M < 4 fall  into  the  domain  of 
large  impactora,  while  cases  where  M > 4 are  in  the  domain  of  small  Impactors. 
The  use  of  the  design  curves  in  each  of  these  domains  will  be  discussed 
separately  here. 


a.  Large  Impactors 

The  procedures  for  the  clamped  and  simply  supported  plates  are  the 
same  but  different  curves  are  employed.  For  simply  supported  plates,  read 
g^OO  from  the  lower  curve  of  Fig,  4 and  e corresponding  to  the  value  of  M 
from  Fig.  26}  for  clamped  plates,  read  g2(n)  from  the  upper  curve  of  Fig.  4 
and  T from  Fig.  27.  Then,  the  maximum  strain  may  be  computed  from  the  formula 


x.max 


■ 8t(n)  e 


vh 


iET 

% ph 


(105) 


where  v • impact  velocity 

h • plate  thickness 
p ■ plate  density 
- plate  flexural  rigidity 
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26.  Selected  Design  Curve  for  Simply 
Supported  Plates  Impacted  by  Large 


Plate  Impacted  by  Large  Iatpactors. 
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Figure  28.  Selected  Deelgn  Curve  for  Impect 
of  Simply  Supported  Rectenguler 
PUte  by  • Smell  Impector. 
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This  procedure  will  be  illustrated  by  an  example. 


Example  #3 


For  a rectangular  simply  supported  plate  with  the  following  properties 


“ 680.  N-a 

°22  " 2490 * K_n 
a * 170  am 

b - 350  am 


h ■ 7 tn 

3 

p - 1.7  g/cm 
E - 8.96  Gpa 
v - 0.3 


impacted  by  a projectile  with  the  following  properties 


R ■ 25  am 


*2  - 1.42  kg 
B - 2.1xl02  Gpa 


v - 0.33 
v » 3.0  m/s 


we  first  compute  the  values  of  the  mass  ratio  and  aspect-orthotropy  ratio 


M - 0.5 


n - 0.451 


From  Fig.  4,  we  find  g^ (0.451)  ■ 0.26,  and  from  Fig.  26,  “ - 0.91  corresponds 
to  M ■ 0.5.  The  maximum  strain  in  the  x-directlon  is 


e,  • (0.26) (0.91) 

XfiBV 


LMKo^oot) 

680  W2 


6.57  x 10 


b.  Gmail  Espactora 

The  design  curve  for  Impact  of  plates  by  small  iapactors  can  handle 

only  simply  supported  boundary  conditions.  To  use  the  design  curve,  read  the 

* 

value  of  generalised  strain  ex  from  Fig.  26  , and  tha  value  of  g(8pin)  from 
Fig.  22.  (Note  that  g(8^,n)  does  not  depend  significantly  on  0p  for  0^  > 7). 
Then,  the  aaximum  strain  nay  be  computed  from  the  formula 


Sou  ” cl  4(1'068)  8(Vn,ep 


(106) 
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where 


6p.  0.09897 


This  procedure  is  illustrated  by  an  example. 

Example  #4 

Consider  the  same  impact  situation  as  in  Example  #3  but  change  the 

impactor  mass  to  m^  * 0.035  kg.  Nov  we  have  14*20 , a small  impactor  case. 

From  Fig.  28,  we  read  e*  ■ 0.025;  and  from  Fig.  22,  g(8  ,0.451)  * 0.35. 

x p 

Computing  6p,  we  get 


8 - (0.09897) 


CO. 

1+0. 


*1221  J. 

0.451  U 

. r {hpitt*.' 

L w- 


(1700X0.007) 

680 


g 2 1/5 

035) 2 


• 6.54 

Therefore,  the  maximum  strain  is 

e - 4(0.025) (1.068) (0.35) (6. 54) 

*•“**  At  680 


1 Pou 

f (1706)T6T007T 


- 6. 79  x 10 


-4 


B.  Construction  of  New  Design  Curves  by  Experiments 

For  beams  end  plates  with  boundary  conditions  different  from  those  con-1 
eidered  In  this  report,  design  curves  may  also  be  constructed  based  on  impact 
experiments.  For  these  structures,  the  parameters  governing  the  impact  re- 
sponse have  the  same  form,  but  the  relationships  (and  thus  the  design  curves) 
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are  different.  For  still  other  structures  (e.g.,  shells,  rods),  new  param- 
eters nay  have  to  be  defined. 


1.  Impact  Strain  Curves 

To  construct  new  design  curves  for  large  inpactors,  one  nay  perform  a 
series  of  intact  experiments  and  measure  the  maximum  bending  strain  which 
occurs  in  the  structure  during  each  impact.  Then,  compute  the  value  of  the 
generalized  strain  according  to  the  appropriate  formula: 


Cb  k 

"max  hv 

c k 


e m £ " 

x y, max  hv 


for  beams 
for  plates 


Plot  one  point  for  each  impact  in  e vs,  M coordinates.  Logarithmic  scales 
are  most  convenient.  Finally,  estimate  a curve  to  fit  these  points.  This 
design  curve  may  then  be  used  to  predict  the  response  of  similar  structures 
having  different  dimensions  subject  to  various  impact  conditions,  as  described 
In  Section  A above. 


2.  Critical  Impact  Velocity  Curves 

In  simple  impact- to- failure  experiments  when  the  strain  Is  not  simul- 
taneously measured,  all  that  is  determined  is  the  lowest  impact  velocity  (or 
critical  velocity)  at  which  the  structure  fails.  However,  this  Is  eufflcient 
to  construct  a design  curve  useful  in  predicting  failure#  due  to  other  impact 
situations. 

For  each  experiment,  the  dlmenslonleee  critical  velocity  v,  defined 
for  plates  in  Eq.  (33)  and  for  beams  in  Eq,  (62),  is  plotted  vs.  the  mass 
ratio  M.  Than  these  points  may  be  connected  to  form  the  design  curve. 
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Mote  that  for  different  impact  conditions,  different  mechanisms  (bending, 
shear,  torsion,  etc.)  may  be  responsible  for  failure;  consequently,  several 
design  curves  may  be  required,  especially  to  completely  describe  a series  of 
widely  varying  impact  experiments. 

Finally,  note  that  since  the  contact  effects  depend  on  the  impact  velocity, 
through  a non-linear  relation,  this  approach  is  not  useful  in  the  small  im- 
pactor  domain. 

C.  Construction  of  Mew  Design  Curves  by  Analytical  Tools 

Using  analytical  methods,  design  curves  for  some  structures  which  are 
not  treated  in  this  report  may  be  constructed.  We  can  offer  a few  suggestions 
on  approaches  to  follow  in  order  to  construct  such  curves. 

1.  Design  Curves  Based  on  Generalised  One-PeRree-of -Freedom  Model 

For  impacts  of  other  structures  In  the  large- impact or  regime,  the  gener- 
alised oae-degreo-of-freedom  model  presented  in  Section  II  may  be  applied. 

To  follow  this  approach,  the  equivalent  structural  stiffness  and  the  strain- 
displacement  factor  d12  must  be  found;  these  may  be  either  derived  from  the 
exact  solution  of  the  corresponding  static  problem  or  computed  numerically 
uaing  soma  approximate  solution  (such  ss  ths  finite-difference  or  finite- 
element  methoda).  Than  the  equation  of  the  daalgn  curve  may  be  derived 
directly  from  Eq.  (8).  For  bsams,  the  form  cf  the  generalized  strain  will  be 
the  earns  ss  In  Eq.  (12);  for  plates,  the  generalised  strain  defined  in  Eq.  (21) 
can  be  employed . For  other  structures,  a new  fora  may  have  to  be  defined. 
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2.  Design  Curves  Baaed  on  Inpact  Calculations 

Several  dynamic  solution  methods  are  available  for  determining  impact 
response  of  structures  on  a problem-by-problem  basis.  These  include  the 
finite-element  method  (in  particular,  NAS TRAN;  see  Ref.  [1]),  Timoshenko- 
type  solution  methods  (application  of  this  approach  to  several  beam  and  plate 
structures  is  outlined  in  [5]),  and  possibly  some  approximate  solutions  based 
on  the  Timoshenko  method  (similar  to  those  presented  in  this  report  for  small 
impactors).  Design  curves  may  be  constructed  by  plotting  the  results  of  a 
few  calculations  using  one  of  these  methods  in  coordinates  of  the  appropriate 
generalized  strain  (e  for  large  impactors,  e*  for  small  impactors)  vs.  the 
mass  ratio  H over  the  range  of  M of  interest.  These  plotted  points  nay  then 
be  connected  to  form  the  design  curve. 
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V.  SUMMARY 

In  this  report,  a method  for  generating  a design  curve  which  predicts 
the  peak  response  of  a structure  subjected  to  low-velocity  impact  has  been 
presented.  Only  a few  experiments  or  analytical  calculations  are  required 
to  construct  such  a curve  for  a given  type  of  structure.  The  importance  of 
such  a curve  is  that  all  impact  cases,  involving  various  impact  velocities, 
structural  dimensions,  and  material  properties,  fall  on  the  same  curve  within 
a variation  tolerable  to  the  designer. 

To  facilitate  our  study,  we  have  divided  all  impacts  into  two  regimes,, 
large  and  small  lmpactors,  each  of  which  exhibits  different  kinematic  behavior. 
Specifically,  oasalve  impactors  have  been  observed  to  produce  multiple  impacts, 
the  peak  response  of  the  structure  depending  on  the  initial  kinetic  energy  of 
the  impactor;  in  contrast,  small  impactors  strike  the  structure  only  once, 
the  maximum  response  being  governed  chiefly  by  the  impactor's  initial  momen- 
tum and  by  contact  effects  at  the  impact  point.  These  differences  in  impact  .. 
response  are  accounted  for  In  the  distinct  forms  of  th«v  design  curves  for  each 
domain. 

Design  curves  have  been  developed  here  for  predicting  the  impact  response 
of  simply  supported  beams  and  simply  supported  and  c lamped-edge  anisotropic 
pistes.  These  curves  have  been  constructed  using  both  results  of  analytical 
solutions  and  data  from  Impact  experiments.  In  addition,  detailed  procedures 
have  been  described  for  using  the  design  curves  presented  in  this  report,  as 
well  as  for  generating  new  curves  for  structures  not  treated  here. 
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APPENDIX  A 


STATIC  SOLUTION  OF  A 
SIMPLY  SUPPORTED  ORTKOTROPIC  PLATE 


In  this  appendix,  we  will  present  the  details  of  the  de- 
rivation of  the  constants  and  d^j  for  the  central  transverse 
impact  of  a simply  supported  orthotropic  plate.  As  explained  in 
the  twin  text,  these  constants  are  needed  for  developing  the  design 
curve . 


I 


The  static  deflection  of  such  a pla'e  due  to  an  arbitrary 
load  distribution  p (x,y)  is,  according  to  Timoshenko  and  Woinoveky- 
Krieger  ^13]  , 


A»* Sm  st  S/f1  n? 


/**#l7r* 


* 


9 


<A1) 


whore 


and 


//  — VjtDk  +2 

r~- 9 


For  a centrally  applied  point  load  P, 


r>  4F*  BE  <■;„  »2L 

p*>»  s,h  2 s,y>  IT 


87 


NADC7C259  -60 


The  maximum  deflection  occuring  at  the  center  of  the  plate  is, 


topis' 


r 


±P*i 

7 T+Oc 


* / 


where 


/*  _ , 4"  r'j  2 2./#  ^ ~ ^ 

~ to  + Jto  » J 


'fe  will  assume  that 

h*  YHty 


(A2) 


so  that 


Therefore,  recalling  that 


we  have 


where 


K1  " P/wl 


ffy)-Z  Z r 


(A3) 


(A4) 


(AS) 


(A6) 


88 


NADC-78259-60 


NADC-7S259-60 

case  of  an  isotropic  plate.  The  advantage  of  this  assumption  is  to 
reduce  the  parameters  for  describing  the  geometry  and  anisotropy  of 
the  plate  to  a single  quantity,  q . For  the  composite  plate  on  which 
the  impact  experiments  were  performed,  the  computed  values  are 
H = D12  + 2 D66  » 1484.6  N-m 
*foxDy  “ *^1^22"  * 1301.2  N-m 
a difference  of  14%. 
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APPENDIX  B 

STATIC  DEFLECTION  OF  A CLAMPED  RECTANGULAR 
ORTHOTROPIC  PLATE  DUE  TO  A CENTRAL  POINT  LOAD 

In  constructing  an  impact  design  curve  for  a particular 
structure  according  to  the  generalized  one-degree-of-freedom  model 
presented  in  the  main  text,  it  is  necessary  to  know  the  relation- 
ships between  load,  deflection,  and  strain  in  the  parallel  static 
problem  for  the  structure  of  interest.  However,  in  the  case  of  a 
clamped  rectangular  orchotropic  plate,  a review  of  recent  literature 
indicates  that  an  exact  solution  for  the  static  deflection  due  to 
a centrally  applied  concentrated  load  does  not  exist.  In  fact,  the 
most  advanced  related  solutions  found  in  the  literature  include  a 
single-term  Ritz  solution  cf  the  case  with  uniform  load  (see 
Lekhnltskli  or  Ashton  and  Whitney  (»0  ) and  an  exact  so- 

lution of  the  simply  supported  case  wit*.  , t :uad  (see  Ambartsumyan 
IvQ  or  Advanced  Composites  Design  Guide  t«0  )• 

Therefore,  in  this  appendix,  the  deflection  and  bending 
strain  of  a centrally  loaded  c lanpssd-edge  rectaugular  orthotropic 
plate  arc  derived  from  the  solution  for  a similarly  loaded  and 
supported  isotropic  plate  due  to  Young  . This  is  accomplished 

by  reducing  the  governing  equations  for  both  problems  to  the  same 
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dimensionless  form  by  making  a simplifying  assumption  suggested  by 
Timoshenko  and  Woinowsky-Krieger  (n]  , and  then  transforming  the 
isotropic  solution  into  terms  of  the  dimensionless  variables  corres- 
ponding to  the  orthotropic  case. 


Young’s  Solution 

In  Young's  solution  for  the  static  deflection  of  a clamped 
isotropic  plate,  the  solution  for  a simply  supported  plate  with  a 
central  concentrated  load  is  combined  with  that  for  a simply  supported 
plate  with  distributed  bending  moments  along  the  edges.  The  edge 
moments  are  then  chosen  so  that  the  deflection  slope  vanishes  at  the 
boundaries.  The  superimposed  solution  consists  of  three  parts, 


w ■ w,  + Wo  + w. 


'1  ’ w2  ' ”3 

where  the  first  term  is  the  solution  for  a simply  supported  rectan- 
gular plate. 


- sM 


where 


- 


z<* 

toira 


The  other  two  terms  are  the  deflections  of  a simply  supported  plate 
with  moments  applied  along  the  pairs  of  dtfgas  y • + b/2  and  x - + a/2, 
respectively. 


/ 
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P«*  ^ „ f-tf£ 

"*  s - -r-gp  <£  /?,„  -■ 

2tt  D *»  fa^ccslu/fr 

[¥**¥  — <*/»  &sA  J 

(s 

and 


tn/rx 

A 


„ P£  _ 

J ' ~ 2 tt*d  tr 


, a»-/ 

» f-d  * 

11  1 


/» 


*tf“7 


't  rrt 

!»£x_^  fmifa  cesk  ^T~] 

the  edge  moments  corresponding  to  and  w-j  are 


W 


jnrrX 


MLu/2  - PS  f-t)  *■/}»  &s  ~ 

* 


(B5) 


itt  which  the  coefficients  and  BQ  are  determined  from  the  condition 
that  the  slope  at  the  boundaries  is  tero.  Young  has  computed  approxi- 
mate values  of  the  first  few  of  these  coefficients  for  several  values 
of  the  aspect  ratio  (b/a). 
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Modified  Young's  Solution 

In  using  the  solution  for  static  deflection  to  construct  a 
design  curve  for  impact  response,  we  are  interested  in  the  relation- 
ship between  strain  and  deflection.  In  Young's  solution  as  described 
above,  differentiation  of  the  term  w^  (twice)  to  obtain  an  expression 
for  bending  strain  leads  to  an  infinite  series  which  is  divergent. 

To  avoid  this  difficulty,  an  alternate  solution,  Navier's  solution, 
which  yields  a bounded  value  of  bending  strain,  is  substituted  for  wj. 
Thus,  we  instead  let 


w,  * 


4 Pa  l 

ir*D 


Transformation  to  Orthotropic  Plate 

The  governing  differential  equation  for  the  deflection  w of 
an  isotropic  plate  subjected  to  a lateral  distributed  load  q(x,y)  is 

^?XZ2/Z  * D 

(*a) 

Letting  x • */a,  y - y/b,  and  v - w/a,  this  equation  may  be  expressed 


l 


'ef'w 


+ 2 


/a\zD4t?  . 9 V 

( V ( v a*  * 
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Likewise,  the  governing  equation  for  an  orthotropic  plate 


n -h  2 H — r *Py 


3V  _ 


(bio) 


n^y  be  expressed  in  terms  of  the  same  dimensionless  variables  as 

- 14W  t = Hi 


* ?v 


(Bll) 


where 


r n 


z-Wffy- 


(M2) 


and  where  it  is  assumed  that  H **  »^xhy  , 

Note  the  similarity  between  eqs.  C®9)  and  (Bll).  The  so- 
lution of  Eq.  (B9)  is  of  the  form 


w - f(  U/hJ  2.  q«J/Dx.  *»y> 


(B13) 


and  the  solution  of  eq.  (Bll)  is  the  same  function  but  involving 


different  dimensionless  parameters. 


w * f(  n.  qa3/i>x.  **y> 


(Bl«) 


Now  eq.  (gl3)  may  be  obtained  by  simply  writing  the  known  solution, 
eqs.  (SI),  (S3),  <SM,  aod  <B?>*  of  thc  isotroPlc  cas*  10  tetta*  of 

the  dimensionless  variables  x,  y and  w.  Then,  substitution  of  n for 
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2 3 ^ 

(a/b)  and  of  qa  /D^  for  qaJ/D  yields  eq.  (814),  which  is  the  solution 
of  eq.  (311).  Finally,  by  returning  eq.  (R14)  so  obtained  to  dimen- 
sional variables  x,  y,  and  w,  we  obtain  the  solution  of  the  orthotrcpic 
equation  (BIO). 

The  terms  to  the  right  of  the  equals  signs  in  eqs.  (B9)  and 
(811)  represent  dimensionless  forcing  functions  for  an  arbitrary  dis- 
tributed load  q(x,y),  and  have  the  form  q ■ qa^/D.  In  the  case  of  a 

*1 

uniformly  distributed  load,  q(x,y)  • qQ,  we  would  have  q - qQa  /D,  or, 
in  terms  of  total  load  P • q0ab,  q • Pa2/bD.  This  suggests  that  in 
the  case  of  a concentrated  load  P the  dimensionless  forcing  function  is 
also  q « Pa  /bD  for  the  Isotropic  case,  or  q • Pa  /bDx  for  the  ortho- 
tropic case.  Therefore,  in  obtaining  eq,  (814)  from  eq.  (013),  Pa2/bDx 
is  substituted  for  Pa2/’oD. 

For  example,  equation  (07)  may  be  written  in  terms  of  dioen- 
tionless  variables  as 


/ jjj7  » 

• a 


rr+bD 


XS2 


60S 


(815) 


/ 


This  represents  a portion  of  the  solution  to  eq.  (013). 

Substitute  q fo?  (a/b)2  and  (Pa2/bOx)  for  (p*2/bD)  to  obtain 

— 4 Pa 2 ^j-4  ^ ftsfifir x c&shtt/ 

/ = fa*  + 

CS16) 

• • 

which  represents  a protion  of  cq.  (814).  Rewriting  this  equation  in 
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terms  of  dimensional  variables  yields  the  final  solution  for  the  ortho- 
tropic case. 


yo  3 ^ tnirx  _ wry 

W, . ^rhr  z Z -iZT-ffTt 

ri&  * * 


(817) 


Similarly,  from  eqs.  103)  and  (B&),  we  obtain 


. -/P*9  _/_  2 
Z Oyr&L  Al  *i/u  *** 


CCS 


torr* 


<ai8> 


/*-/ 


to'  c ^ v- - •wer*  jo 

3 SrrHO*  fjf  5 ^ 


■*  2tt*i,d„  uj  tr  »***<*,  * 

/t^  ^ 7^  -*»&**«* *«*  7%f] M a, 

^ ^ ? 22  y~/)^  y.  mn 

J wo,  w?b’£^*u- 

By  assembling  cqs.  (El?),  (818)  and  (8191  according  to 
eq.  (81),  we  obtain  the  solution  for  the  deflection  of  a rectangular 
clamped- edge  orthbtropic  plate  subjected  to  a centrally  applied  con- 
centrated load.  The  coefficients  Apj  and  B^,  which  are  functions  of 
(a /b)  in  cite  isotropic  case,  ere  now  functions  of 

Evaluating  the  deflection  solution  at  the  center  of  the 
plate  (x  • o,  y » 0),  v«  obtain  the  maximum  deflection  of  the  plate 


which  may  be  expressed  as 


W - I P 


: W p 


(820) 
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SO  that 


IS  - jk.P* 
K/  *3 


A 


(B21) 


wher«? 


w- 


2tt 


Z [Am  ^ ** 


’Vf  » 

+d» 

toSh  pv»  J 


(B22) 


and  where  we  now  have 


««* 


M7T 


(B23) 


The  maximum  bending  strain  occurs  at  the  point  directly  opposite  the 
load  and  may  be  expressed  as 


E* 


h a 

'***-  b 0K 


Wp 


(B2A) 


-2. 


&>sA  s* 


ys.  B*__  7 

f * £>st/dh,  J 
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APPENDIX  C 

COMPUTER  PROGRAM  FOR  CALCULATING 
TIMOSHENKO  SOLUTION  OF  PLATE  IMPACT 


This  appendix  presents  a listing  of  the  FORTRAN  computer  program, 
•named  SMINC3,  for  calculating  the  Timoshenko  small-increment  solution 
of  the  central  transverse  impact  of  a rectangular  orthotropic  plate, 
Iq,  ( 40). 

The  program  as  listed  is  set  up  to  solve  an  impact  problem  having 


the  following  parameters: 

v ■ 2.45  m/sec 

a x b x h - 171  x 349  x 7.1  on 
n^  - 0.660  kg 

Oj  » 0.816  kg 

The  plate  bending  stiffnesses,  computed  according  to-  the  Whitney- 
Pagano  (1970)  anisotropic  plate  theory,  are: 


°11  " 680,4  N”Q 
D12  - 473.2  N-u 

D22  • 2488.8  N-« 

" 505.7  N-» 

°16  * °26  “ **49,1  N"° 

A44  " A55  " 2,687  * 10?  N/o 


The  plate  is  treated  as  specially  orthotropic  in  banding  and  the 
and  D2g  terms  are  neglected. 

The  program  Itself  calculates  the  value  of  the  Hertslan  contact 


stiffness  constant  kj  bated  on  supplied  values  of  isotropic  elastic 


I 


100 


N ADC -78259 -60 

properties  of  the  plate  and  iopactor.  The  values  used  are: 


plate  (Gr/Ep) 


impact or  (steel) 


f E - 9.8  : 
- 0.3 

f E"  2.1 

- 0.3 


x 109  N/m2 


x 1011  N/o2 


The  contact  radius  of  the  impactor  is  25.4  tmn. 
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C SMINC3  — PLATE  IMPACT 

COMPUTATION  OF  TIMOSHENKO-TYPE  SMALL-INCREMENT  SOLUTION  OF  HERTZIAN 
C IMPACT  OF  A SIMPLY  SUPPORTED  ORTHOTROPIC  RECTANGULAR  PLATE 
C MODIFIED  SUN  fi  CHATTOPADHYAY  METHOD 

C DT  - TIME  INCREMENT 

C NOT  - NO,  OF  TIME  INCREMENTS  CALCULATED 
C TOL  - TOLERANCE  OF  ERROR  IN  NONLINEAR  SOLUTION 
C V - IMPACT  VELOCITY 

C XA»  YB,  H - DIMENSIONS  OF  PLATE  IN  X,  Y*  Z DIRECTIONS 
C Ml,  M2  - MASS  OF  PLATE,  IMPACTOR 

C El,  E2  - ELASTIC  MODULUS  OF  PLATE,  IMPACTOR 

C PRl,  PR2  - POISSON'S  RATIO  OF  PLATE,  IMPACTOR 

C NHM,  NHM  - NUMBER  OF  VIBRATION  MODES  CONSIDERED 

C K2  - HERTZIAN  CONTACT  STIFFNESS 

C XNX,  YNY  - INITIAL  STRESSES  IN  X,Y  DIRECTIONS 
REAL  L,MI,M2,K2,L11,L12,L13,L22,L23,L33 
DIMENSION  FI5000) «GW(5000) ,GXX( 5000! ,GYY 15000) 

£**«»****  PROBLEM  DATA  ********************************************* 
NHM*25 
NHN*25 
NOT- 100 
TOL-l.OE-05 
DT*1.0E-06 
XNX*0. 

YNY*0. 

011*680.4 

012*473.2 

D22*2488.6 

066*505.7 

A 44*2,6876  07 

A55*A44 

XNX*0. 

YNY»0. 

XA*0, 171 

YB-0.349 

H*0.0071 

Nl-0.660 

M2*0.816 

Rl*0.0254 

PRl*0.3 


PR2-0.33 
£1*9.86  09 
£2*2. 10E  ll 
V*2.45 


£****•**••••*#**•*******«#«*•***♦***•**********************«********* 


Z-H/2. 

PI  *3. 14159265398 

1001  FORMAT l4X,'N',aX,'TINE*,l4X,«fORCE*,12X,' APPROACH* ,10X, 'DEFLECTION 
l',10X* 'STRAIN  X',10X, 'STRAIN  Y' ,10X, 'STRAIN  XY' ) 

0EL1*I l.— PRl**21/6l/PI 
0EL2-1 l *-PR2*42l /62/P I 
K2-4.*SQRT(Rn/(D€LU0EL2)/PI/3. 

CAY-PI**2/12. 

A44*A446CAY 

A55-A554CAV 

B*0TM2/3.O/N2 


C-  4.0/Ml 
P-M1/XA/YB 
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PRINT1001 

CALCULATE  NATURAL  FREQUENCIES  OF  PLATE 
D090K=1,N0T 
GWCK)*0. 

GXXCKJ-O. 

90  GYY(X)=0. 

D0247M=i»NHM,2 

TERMA*M*Pl/XA 

TERMA2*TERMA**2 

L13=A55*TERMA 

D0247NS1 »NHN,2 

TERM8=N*PI/YB 

TERM82=TERM8**2 

LU*Dll*TERMA2*  D66*TERMB2*A55 
L12“( D12+D66)*TERMA*TERMB 
L22»066*TERMA2+D22*TERMB2+A44 
L23=A44*TERMB 

L33=CA55*XNXI*TERMA2M  A44*YNY  WERHB2 
Q=L11*L22-L 12**2 

DETERM=Q*L33+2.*L12*L23*L13-L22*L13**2-L11*L23**2 

AC=t  L12*L23-L22*L 13) /Q*TERMA 

BC=(L12*L13-Lll*L23)/Q*TERMB 

0MEGA2=DETERM/Q/P 

OMEGA= SORT  COM EGA2) 

IF  ( CM. EQ.U.AND.lN.EQ.l)) PRINT  1002, OMEGA 

1002  FORMAT ( # OMEGA, 1,1  = »,E15.8) 

TERM1»DT*0MEGA 

C2-1.0 

D0247K*1 * NDT 
C1*C2 

C2=C0SCK*TERM1> 

AAA=CCl-C2)/0MEGA2 

CONSTRUCT  TABLES  OF  SUMMATION  FUNCTIONS 
GWCK)*GWCK)*AAA 
GXXIK)*GXX(K) ♦ AAA* AC 
247  GYYCK)«GYYCKJ*AAA*BC 
P»B*-C*GWU) 

Q«l./3. 

S-2./3. 

D0100N*1,NDT 

CREATE  HEADINGS  FOR  DATA  AT  TOP  OF  EVERY  PAGE 
CF(M0DCN,6U.NE.0)G0TQ70 

1003  FORMAT ( 1H1 ) 

PRINT1003 

PRINT1001 

70  T»N*DT 
A*V*T 

IFCN.EQ.UG0T014 

8UM-0. 

NMi-N-l 

SUM1-0. 

SUM-0. 

005J-I.NMI 

K-N-J 

R*K 

SUM  I -FU) -SUM  I 
SUN-SUM*R*SUMI 
5 BUM«BUM*F(J)*GW<K*1) 
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A*A-C*BUM 

A»A-DT**2*(2.*SUM-SUMI/3,I/M2 
COMPUTE  SOLUTION  FOR  FIN)  USING  NEWTON’S  ITERATIVE  METHOO 
G0T013 

14  FNEW=SQRTU**3)/K2 
13  FN=FNEW 

B»SIGNU.0»FN1/K2**2-3.0*A*P**2 
FFN»P**3*FN**34B*FN**2 
1 «-3.0*A**2*P*FN-A**3 

FPFN*3,0*P**3*FN**242.0*B*FN«-3.  0*A**2*P 
82  FNEW*FN-FFN/FPFN 

ERROR*ABS( (FNEW-FNJ/FN) 

I FI ERROR-GT.TOL I G0T013 
IF! N*GT« l JG0T010 
CHECK  THAT  Ftl)  IS  POSITIVE 
IFIFNEW.GUO.IGOTOIO 
FNEW=-FNEW 
G0T013 
10  F(Nl*FNEW 

IFIFNEW,LT.0,)F(N1*0. 

CALCULATE  STRAIN  ANO  DEFLECTION 
W*0. 

PSIXX=0. 

PSIYY-O* 

D021J=l»N 

W*W«F(J1*GWIK) 

PSIXX*P$IXX»F( J)*GXXIK> 

21  PSIYY«P$IYY»FIJ)*GYYIK» 

w*w*c 

EPSX—PSIXX*C*Z 
EPSY— PSIYY*C*Z 
GAMMA*0. 

ALPHA* SIGN! l«0*FNEW|*| ABSIFNfcWJ/K2l**S 
PRINT1000«N»Tf FINlf ALPHA»W»EPSX«EPSY»GAMNA 
100  CONTINUE 
STOP 

1000  FORMAT !lXtI4»7l 3X*E15«8I ) 

FND 
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APPENDIX  D 

SIMPLY  SUPPORTED  BEAM 
CENTRALLY  LOADED  BY  HALF-SINE  PULSE 

In  this  appendix  are  derived  expressions  for  the  deflection 
and  bending  strain  in  a simply  supported  beam  subjected  to  a half- 
sine pulse  laterally  applied  at  midspan. 

The  deflection  of  a simply  supported  beam  due  to  an  arbitrary 
force  F(t)  applied  at  midspan  is 

x tzL 

11  a 

erne 

(Dl) 


„ 2_  & 


W(Xjt)  JOJ&)  —————  Sth 

r^i^xs  &ji  L. 

J Ffr)  s/'/t 


where 


6Jj  = fir^aVt 


We  yUl  consider  the  case  where  the  applied  force  is  repre- 
sented by  a half-sine  pulse  with  respect  to  time. 


Fd)'- 
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Then  the  solution  for  each  time  domain  (during  and  after  the  pulse) 
may  be  written  in  terms  of  the  convolution  integral 

v&jt)  - ^ » °T 

iHjis  L 

t 

f sin  fr  Sin&J;(t-r)dT,  0< 


f Hir 

f sin  ft  sin Ui  (t-rjdrj  -fay  o<t<  r/f 

J6 

d . , 2P  d-/)^  „■  irrn_ 

*r*'rj  fAL  (SXr  Mi  L 

,v? 

Sin  f T sin  tyfy  -r)drj  -6r  £>lr/i 


i: 


Evaluating  the  integral,  the  deflection  during  the  pulse  is  given  by 

Jsfo  it-  siHUjty/A*-  o<t<ir/f 

ft  * i/iJ,  (DS) 

2 

If  g*  i t however,  the  denominator  of  this  expression  is  zero,  and 

2 

the  result  is  undefined.  Xu  this  case,  ^ » i w,  - 0uj  * $;  for  a 
steady-state  vibration  problem,  this  would  correspond  to  a resonant 
forcing  function.  For  this  special  case,  the  following  term  must  be 
substituted.  • . 

opr  * I*.  /~/i~sr  / . 3.. 


Vfi  ~ 


. iirx  /-/)**■ 

V*£l  * 


, ftsm  ft 
2,  / * 

— t cos  ft  +4r  &sftsm2f?) 


<D6) 
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After  the  pulse,  the  deflection  is 


^■iL.  m iF- 

z *~ 


/•/,$, s0  (ft  - i *)  <' 

• Jsm  6J(t  + sm  fr/t-  ir/f)J  ^ t>^f 

2 

unless  6 ■ i , in  which  case  the  ith  term  must  be  replaced  by 

M 

w 2PL3_  s/hlm~  (z£LL  tc»s  ft 

7f4£I  L (D8) 


The  strain  due  to  bending  is 


£=  ~Z 


2JV 

9x*- 


(09) 


The  maximum  value  of  this  strain  is  reached  at  midspan.  If  we  define 
• dimensionless  strain  as 

A rraEI  P 
e = h TL 

then  the  solution  In  tertas  of  this  quantity  is 

r » 


(010) 


jt  ~-(sm  ft  - i O < £<rr/f 

!W  V ' * S (Oil) 

X<  "jCt/*******# [t-r/fj),  t>ir/f 

(D12) 
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2 

unless  8 ■ i , in  which  case  the  ith  term  must  be  replaced  by 


(sin lft- 


(D13) 


\~2 


(DU) 


For  most  problems  involving  small  impactors  (large  M),  the  contact 
duration  is  short  compared  with  the  fundamental  period  of  beam  vi- 
bration. Therefore,  it  is  of  interest  to  determine  the  peak  bending 
strain  when  the  pulse  duration  is  short.  It  is  observed  in  this  case 
that  the  peak  strain  occurs  after  the  pulse,  that  is,  in  the  domain 
of  eq.  (D12),  which  may  be  written  as 

e = Z sh[l  V/*-  f)/z]cos  ^ 

H*. SP'l 

This  quantity  reaches  its  maximum  at 

t 7r 


t » 


when  It  has  the  value 


A1  4 " 2/3 

which  is  a function  of  6 only.  This  function  is  plotted  in  Fig.  Dl. 


£ « 
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Maximum  . .menslosless  sc  rain  eBt33t  vs*  frequency  ratio  B for  a simply  supported  beam 
subjected  to  a bnlf-nirtc  pulse. 
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APPENDIX  E 

DETAILED  RESULTS  OF  IMPACT  EXPERIMENTS 

In  this  appendix,  data  from  all  intact  experiments  conducted  at 
Dr  axel  and  also  some  reported  in  the  literature  are  presented  in  both 
tabular  and  graphical  form*  for  completeness,  several  figures  which 
appear  in  the  text  have  been  included  in  this  appendix* 

U Results  of  Large-Iapactor  Strain  Measurement  Experiments 

Strain  measuring  ispact  tests  have  been  performed  on  both  aluminum 
and  composite  ( graphite- epoxy  ) plates*  A total  of  475  experiments 
were  performed  on  specimens  with  clamped  edges  and  573  on  specimen*  with 
ajjqply  supported  edges* 

Index  to  Large  lapaotor  Design  Curves 

Platea  Simply  Supported  Plates 
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X Impact  Experiments  on  Aluminum  Plates 
One-Degree-of-Freedom  Model 

Figure  B3«  Design  curve  and  experimental  data  for  Impact  of  a 
clamped  aluminum  plate  (1/n  ■ 0*25). 
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Figure  14*  T)c .ign  curve  and  experimental  data  for  Impact  of 
a clamped  aluminum  plate  (n  *9.0). 
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X Experiment* 


Figure  123#  Design  curve  and  experimental  data  for  Impact  of  a 
•Imply  supported  aluminum  plate  (n  • 1.0). 


Mass  Ratio,  M * 

«•»  One-degrec~of-fresdo» 

K Experiment 


Figure  124*  Design  curve  and  experimental  data  for  impact  of  a simply 
supported  plate  (n  * 4.0). 


Mass  Ratio,  M ■ m^/n^ 

■t 

■—  — - One-degree~of -freedom  model 
X Experiment 

Figure  S27*  Design  curve  and  experimental  data  for  impact  of  a 
simply  supported  aluminum  plate  (1/n  * 0.111). 
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Figure  S41. 
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2.  Results  of  Imp  act-to-F allure  Experiments 

Each  specimen  was  Impacted  at  a gradually  increasing  velocity 
until  damage  was  detected  either  by  inspection  or  by  use  of  a hand  held 
ultrasonic  transducer  and  oscilloscope  system. 

The  theoretical  failure  strain,  £f,  used  In  constructing  the  one- 
degree-o f- freedom  model  design  curve  is  the  surface  strain  obtained 
when  the  strain  in  the  outermost  0*  lamina  reaches  the  failure  strain 
supplied  by  the  manufacturer  (0,0112). 


Impact 

Mass, 

o2  (kg> 
0.1095 


0.0566 


MAOC -78259 -60 


TABLE  S5' 

:t  failure  experiments  performed  on 

COMPOSITE  BEAM  SPECIMENS 


Mass 

Ratio, 

Impact 

Velocity, 

v (m/sec] 

Failure 

Occurred 

Dlmensionlcs: 

Impact 

Velocity, 

v * vh/a^ 

0.0562 

5.63 

No 

0.00333 

6.32 

No 

0.00374 

6.77 

Ni 

0.00400 

6.95 

No 

0.00411 

7.90 

No 

0.00467 

7.00 

Yes 

0.00414 

7.35 

Yea 

0.00435 

7.33 

Yes 

0.00445 

7.65 

Yes 

0.00452 

8.25 

Yes 

0.00486 

8.66 

Yes 

0.00512 

' 

G.1G39 

9.31 

No 

0.00350 

9.42 

No 

0.00537 

10.30 

No 

0.0060S 

11.20 

No 

0,00662 

11.13 

Yes 

0.00658 

11.27 

11.47 


Yes 

Yes 


0.00666 

0.00678 
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TABLE  85  (continued) 


Mass 

Ratio, 


Failure 

Occurred 


Dimensionles 

Impact 

Velocity, 

v * vh/a^ 


.0566 


0.0280 


0.01434 


0.1089 


0.220 


0.429 


11.72 

Yes 

0.00693 

12.97 

Yes 

0.00767 

13.03 

No 

0.00770 

15.45 

No 

0.00913 

17.59 

Yes 

19.17 

■an 

0.01133 

20.12 

Yes 

0.01190 

23.07 

Yes 

0.01364 

24.07 

Yes 

0.01423 

24.10 

Yes 

0.01425 

19.05 

No 

0.01126 

21.17 

No 

0.01252 

21.97 

No 

— 

22.66 

Yes 

0.01360 

23.39 

Yes 

0.01383 

26.24 

Yes 

28.42 

Yes 

0.01680 

33.43 

Yes 

0.01976 

NADC-78259-60 


■i 


No  Bending  Failure 


0.02  ** I-.- .. ■> — — — l. -J 

0.02  0.05  0.1  * 0.2  0.5  1.0 

Hass  Ratio,  H 

—— — Cr Ideal  Velocity  fair  Shear  Failure 

Critical  Velocity  for  Bending  Failure 
Experiments  on  Graphite-Epoxy  Beams: 

H Failure 
• Ho  Failure 

Figure'E^5*  Critical  impact  velocities  for  shear  and  bending 

failure  and  experimental  data — a bending-dominated 
case. 
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HH-  jSf 51=2=3:,., 


• Simply  Supported  Plates 


• - *-*-M  •••-’  -«*•*.*.  •«*«»•»  <» 

U*  *«»•  **•»  .*  i>«<t  »>•>•««  ».« 


1.5  /L  2.5  ,3  / £ J Ji  9 


Ifi  lie  i 2.5  3 4 5 i 7 8 

Impact  Velocity  Curve  for  M 


Figure's^,  Critical  Impact  Velocity  Curve  for 
Clamped  and  Simply  Supported  Plates 
Impacted  by  Large  Masses. 
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3.  Results  of  Small  Inpactor  Strain  Measuring  Experiments 

Plate  specimens  were  Impacted  by  l"(25na)  radii  aluminum  projectiles 
fired  from  an  air  gun.  Impactor  velocities  were  recorded  using  a photo- 
diode system.  Plate  strain  was  recorded  using  type  EA-13-062TT  rosettes 
and  a type  565  dual  beam  oscilloscope. 
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TABLE  S7 

DATA  FROM  EXPERIMENTS  ON  SIMPLY  SUPPORTED 


BEAMS  IMPACTED  BY  SMALL  MASSES 


Beam  Mass, 

Dimensions, 

and  Material 

Mass 

Ratio, 

M*m1/m2 

Impact 

Velocity, 

v (m/sec) 

Contact 

Stiffness, 

k2(N/m3/2) 

Generalised 

Strain, 

e ■ ea2/hv 

e 

/T  ' 

0.389  kg 

192  x 16.x  16  mm 

steel 

27.0 

6.37 

2:466xl010 

0.0795 

0.0184 

II 

13.17 

ft 

0.0814 

0.0175 

13.8 

It 

0.1229 

0.0315 

It 

12.19 

It 

0.1489 

0,0369 

6.85 

5.49 

ft 

0.2388 

0,0739 

It 

8. 53 

It 

0.2512 

0.0743 

3,54 

6.28 

It 

0.3224 

0.1123 

M 

9.33 

II 

0.3493 

0.1171 

0.601  kg 

305  x 16  x 16  oa 

steel 

41.8 

10.88 

M 

0.0800 

0.01142 

II 

14.87 

II 

0.0790 

0.01092 

23.8 

7.92 

It 

0.1033 

0.01703 

M 

11,06 

It 

0.1076 

0.01716 

11.2 

7.50 

♦I 

0.1766 

0.0341 

M 

12.37 

ft 

0.1800 

0.0330 

5.83 

6.16 

1 J 

0.2478 

0.05S4 

It 

8.11 

It 

0.2662 

0.0579 

0.1095  kg,  aluoinuc 

1.93 

4.39 

II 

0.475 

0.1689 

152  x 20  x 12.7  rea 

H 

: 

6.92 

II 

0.462 

0.1569 

203 
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TABLES?  (continued) 


Beam  Mass, 
Dimensions, 
and  Material 

Mass 

Ratio, 

K^m^/m^ 

Impact 

Velocity, 

v (m/sec) 

Contact 

Stiffness 

k2(N/m3/2) 

Generalized 

Strain, 

e ■ ca2/hv 

C 

/b* 

0.601  kg 

305  x 16  x 16  mo 

steel 

41.8 

10.58 

1.744xl010 

0.0760 

0.01166 

n 

12.71 

II 

0.0761 

0.01146 

23.8 

9,42 

It 

0.1242 

■89 

II 

12.53 

II 

0.1235 

0.0208 

11.2 

5.33 

It 

0.1534 

0.0328 

It 

12.01 

It 

0.0317 

5.88 

5.39 

II 

0.2345 

0.0570 

H 

7.71 

If 

0.2799 

0.0656 

0.360  kg 

185  x 16  x 16  ma 

steel 

25.0 

11.31 

2.466xl010 

0.0368 

0.0203 

It 

14.42 

It 

0.0978 

0.0223 

12.8 

8.11 

t» 

0,1523 

0.0421 

It 

10,70 

It 

0.0442 

6.34 

8,20 

It 

0.2563 

0,0814 

It 

10.52 

H 

0,2434 

0.0754 

3.27 

5,73 

*1 

0.3252 

0.1221 

It 

| 

tt 

0.3203 

0.1155 

0.1074  kg,  steel 

98  x 15  x 9 am 

IH9E9I 

t« 

0.458 

n 

■Hi 

6.40 

It 

0.465 

0,2147 
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Bean  Mass, 

Mass 

Dlnenslons, 

Ratio, 

and  Material 

M-n^/Oj 

3.66  kg 

6.44 

762  x 25.4  x 25.4 

am 

•tcel 


56.9  kg 

2900  x 50  x SO  oa 
steel 


3.66  kg,  steel 
724x25.4x25.4  net. 


11.41 

6.34 

55.0 


TABLE  SB  . 

DATA  FROM  GOLDSMITH  (1960) 


Inpact  Contact  Generalized 
Velocity,  Stiffness  Strain, 

v (a/sec)  kjOl/n^*)  t • ea^/hv 


0.405 

1.338  ” 

2.54 

0.408  '» 

1.271  " 

1.777  " 

2.40  4.894X101 

2.47 

UWMMPMMM  MMMMIMHHMM 

3.00  " 


2.40 

2.40 

1.725 


2.170  kg,  steel 
762x19x19  am 

266. 

45.7 

1.447  kg,  steel 
762x19x12.7  tan 

177.2 

45.7 

0.965  kg,  steel 
762x12,7x12.7  mm 

118.2 

45.7 

0.482  kg,  steel 

59.1 

45.7 

1. 744x10* 


1.233x10 


Lo 

.286 

.311 

0 

.520 

.564 

0 

.584 

0. 

148 

.247 

0 

.237 

0 

.260 

0 

,239 

■ 

.365 

p 

.420 

■ 

.0802 

0 

.0268 

■ 

.0359 

G 

.0713 

.0563 


0.0500 


0.0510 


0.1278 


0.1237 


0.1239 


0.00932 


0.01855 


0.0)740 


0.01886 


0.01676 


0.03173 


0.04107 


0.01065 


0.001486 


.1578  0.00 
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M 

* /0 

ATj  4K/s 

? 

** 

Ft 

/L  9 

<9000 

2i,B 

9.38 

0.097/ 

3400 

K.t 

47.38 

0.0935 

/A? 

3000 

Z9.B 

3.72- 

6,  O'/S 

/&60 

/T.7 

3.50 

o.ossf 

k_.  

/too 

77.1 

3,  SO 

0.029? 

gooo 

U>8 

3.00 

a.468 

3600 

19.7 

3.75 

OtO&9 

FI 

323 

3000 

26,8 

46$ 

0,6/3! 

4/000 

26.8 

6,029! 

* 4 

i 

. „ _ 

/too 

!$,3 

5.5? 

0,0/03 

<?& 

4 

18.3 

W^IMWEEM 

Data 


M 


*f.67 


/SS 


{70 
36 * 
63/ 
730 
lot 
7/60 
1263 
{*82 
332 

03 

IT4/ 

973_ 

/Q/7 

/tfi 
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TABLE  E10 
Fr:^  Schwieger 


4T  w/*  $ 6 

0,26  /S,7  0.0/28 


0,50 

z/,s 

0,0388 

0,75 

23,3 

0.0358 

AM 

25,0 

0,63*9 

/.iS 

25,8 

0.033 / 

/,50 

Z&.7 

0.0321 

/.75 

21.  b 

0.03/3 

2,00 

28.3 

0.0308 

6,60 

3i,S 

0,0233 

6, IS 

3&5 

/.oo 

36,(7 

0,0213 

/.25 

38.1 

0,  0/90 

/.So 

39,7 

0.0/92 

ATS 

W,f 

0.0/61 
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Figure  B52#  Strain  Curve  for  Impact  of  Simply  Supported 
Rectangular  Plate  by  a Small  Impactor. 
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Scale: 
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a)  Simply  Supported 
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B/io/ysiH 
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Scale: 
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b)  Clamped 


Plgura  K53*  Impact  of  Specimen  B3,  M - 0.21,  Impact  Velocity  - 1.73  a/a 
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b)  Clamped 


Figure  E54:  Impact  of  Specimen  B3,  M « 0.21,  Impact  Velocity  • 2.45  a/sec 
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a)  Simply  Supported 
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b)  Clamped 


figure  E55:  Impact  of  Specimen 


B7»  M - 2.38*  Impact  Velocity  • 1*73 


N ADC -78259 -60 


1000  II  — /dlv 

Tmn 

vertical 
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Figure  E56:  Impact  of  Specimen  B7,  M - 2.38,  Impact  Velocity  - 2.49  m/sec. 
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Appendix  F 

The  Effect  of  the  Approximation  H » ^n°22  on  t*ie 
Predicted  Impact  Response  of  an  Orthotropic  Plate 

To  simplify  the  analysis  of  the  bending  and  vibration  of  specially  ortho- 
tropic plates,  Timoshenko  and  Wolnowsky-Krleger  [13]  have  suggested  the  approx- 
imation 

H ■ fllD22  <n> 

where  H is  defined  as 

H 5 Di2  + “ee 

An  advantage  of  using  this  assumption  is  that  the  equation  governing  bending 
of  an  orthotropic  plate  may  be  transformed,  through  a simple  change  in  variables, 
into  an  equation  similar  to  the  equation  governing  bending  of  an  isotropic  plate. 
By  reversing  the  transformation,  a solution  for  the  deflection  of  an  isotropic 
plate  may  be  converted  into  the  corresponding  solution  for  an  orthotropic  plate. 
This  is  done  in  Appendix  B for  a centrally  loaded  rectangular  plate  with  clamped 
boundary  conditions;  the  approximation  (FI)  is  instrumental  in  this  instance 
because  no  general  solution  exists  for  the  orthotropic  case. 

Like  all  approximations,  however,  application  of  Eq.  (FI)  introduces  some 
error  into  the  solution.  In  this  Appendix,  we  investigate  the  influence  of  this 
approximation  on  the  strain  predicted  by  the  one-degree-of -freedom  model  of 
impact  of  a simply  supported  rectangular  orthotropic  plate.  This  la  done  by 
deriving  the  relationship  between  the  generalised  strain  and  the  mess  ratio  M, 
both  vlth  and  without  the  use  of  Eq.  (FI),  and  then  comparing  the  results. 
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Th«  static  dsflection  of  a simply  supported  rectangular  orthotropic  plate 

due  to  a central  concentrated  load  P is  [10] 

. gig  j airy  . orr  . nir 
8 in sin  — r*-  sin  — sin  — 

v<*,y> -j — h—s* 

abr  m n m_  D + 2Hs_n_  + D 
a4  D11  + ^ b4  °22 

where  a and  b are  the  plate  dimensions  in  the  x and  y directions.  If  we  define 
n aa  in  Bq.  (16)  and 

r - H 


fv 


22 

than  the  maximum  deflection!  which  occurs  at  the  center  of  the  plate,  may  be 

expressed  as 


- P/K, 


share 


e-  DU  b 


^3  P^n) 


(P2) 


*!<«»*!>  ■[[  f-  ».n 


odd 


n aa 


B_ 


4 2 2 4 o 

a 2m  n nt  + a n 


Similarly,  the  maximum  x-dlrection  bending  strain,  which  also  occurs  at 
the  plate  center,  may  be  expressed  as 


2 Pha 

W *2  “u  * 

2h  a R, 

• 1 

* BUb 


P2U,n) 


P2(C.n)  wmx 


where 


*2U.n)  - ll  f— 

an  an 


m,n  add. 
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Therefore,  the  strain-deflection  constant,  defined  in  Eq.  (6),  is 

w D11  b 1 

d12  “ 2h  a F2(c7nT  (F3) 

To  apply  the  one-degree-of-freedom  impact  model  to  this  case,  we  substitute 
Eqs.  (F2)  and  (P3)  into  Eq.  (8)  with  e • 0 and  define  the  generalised  strain  e 


as  in  Eq.  (21)  to  obtain 


G1(c,n) 


where 

F,(c.n) 

G.(c.n)  - — 

1 ^(S.rO 

Note  that  the  approximation  (Pi)  corresponds  to  setting  c * 1.  Thus,  if 
this  assumption  is  used,  then  the  approximate  relation  between  and  Ht 
corresponding  to  Eq.  (22),  is 


”1 

e.  «■ 


G. (l,n) 


Therefore,  the  correction  required  to  adjust  the  approximate  generalised  strain 
given  by  Eq.  (F5)  to  the  exact  value  given  by  Eq.  (P4)  is 


H 


G (t,n) 

B<ex)  “ G^l.n)  “ 1 

This  error  is  plotted  vs.  t for  a few  values  of  n In  Fig.  FI.  No*a  that 
the  magnitude  of  this  error  does  not  exceed  17X  for  the  rsnge  of  c vslues 

calculated  and  is  laaa  than  9X  for  moat  lay-ups  commonly  used  In  actual  aircraft 
atructurea,  which  lie  In  the  range  0*60  < C < 1.6.  (Lay-ups  having  € valuea 
outalda  this  range  [e.g.,  C • 0.31  for  unidirectional  AS/ 3501-6]  rarely  have 
practical  application  aa  plate  atructurea.)  Further,  for  the  range  of  C valuea 
of  the  plata  specimens  uaed  In  the  impact  experiments,  this  error  Is  always  laaa 
than  6.5X. 
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In  Table  FI,  the  actual  values  of  n»  and  the  error  in  predicted  generalized 

strain  E(e  ) for  the  plate  specimens  are  summarised.  The  largest  value  of  error 
X 

among  all  of  the  plates  is  only  5.4X,  Therefore,  ve  conclude  that  the  approxi- 
mation of  Eq.  (FI)  does  not  contribute  a significant  error  to  the  design  curve 
for  Impact  of  simply  supported  plates  by  large  iapaetors.  He  may  infer  that  a 
similarly  small  correction  is  required  for  orthotropic  plates  with  other  boundary 
conditions . 

Table  FI 


Error  in  Predicted  Generalized  Strain 
for  the  Plate  Impact  Specimens 


Plate 

Humber 

X 

1 

M 

81 

1.430 

0.623 

-3.6X 

B2 

0.164 

-4.2 

B3 

2.39 

-3.7 

84 

0.625 

-3.6 

85 

0.074 

-4.0 

86 

5.29 

-5.4 

87 

1.38 

-3.2 

88 

0.282 

-4.1 

B9 

0.625 

-3.6 

FI 

1.251 

0.641 

-2.2 

F2 

0.168 

-2.6 

F3 

2.45 

-2.3 

F4 

0.641 

—2.2 

Hi 

1.270 

1.150 

-2.1 

H2 

0.302 

-2.7 

H3 

4.41 

-3.3 

H4 

1.150 

-2.1 
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NAVSEASTSCOK,  Uashlt«ton,  D.C.  20352 

(Attn:  Coda  035,  Hr.  C.  Pohlar).  ..............  1 

MAV5EC.  Byattavllle,  HD  20782 

(Attn:  Codo  0101103,  Hr,  W.  Cranes).  • 1 

NAVSHIPRAHOCEN,  Bethcsda,  KD  20034 

(Attn:  Code  173.2,  Mr.  W.  P.  Cauch).  • ••••.«••«••  1 

NAVSBtPtANDCEN,  Annnpolia,  KD  21402 

(Attn:  Coda  26?o,  Mr.  B.  Edelaeein).  ..*...•••••,  1 

NCL,  White  Oak,  M>  20910 

(Attn:  Mr.  P.  R.  Barnet). . . 1 

MU,  Washington,  D.C,  20375 

(Attn:  Dr.  1.  Wolock).  . , I 

OMR,  Washington,  D.C.  20362 

(Attn:  Dr.  N.  Perrons)  I 
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Gcvenment  Activities  (Cont.) 

PLASTEC,  Picatinny  Arsenal,  Dover,  KJ  07801 

(Attn:  Librarian,  Bldg.  176,  SARPA-FR-W-D  and  Mr.  B.  Babiy).  • 2 
Scientific  & Technical  Information  Facility,  College  Bark,  MD 


(Attn:  NASA  Representative).  1 

TJSAAVMATLAB , Fort  Bustla,  VA  23603 

(Attn:  Mt.  R.  Beresford)  . . . , * o • * • • . 1 

SEAMATRESAGj  Watert own,  KA 

(Attn:  Dr.  E.  Lenoe)  1 

0EARES0FC,  Durham,  NC  27701  1 


Kon-Govarnwmt  Agencies 


Avco  Aero  Structures  Division,  Nashville,  TN  37202 
(Attn:  Mr.  W.  Ottanville). 

Bat telle  Columbus  Laboratories,  Metals  and  Ceramics  Information 
C«ut«r,  505  King  Avenue,  OH  43201. 


Bell  Aerospace  Company,  Buffalo,  NY  14240 

(Attn:  Zone  1-85,  Mr.  F.  M.  Anthony)  ..••••••••••• 

Ball  Helicopter  Company,  Fort  Worth,  IX  76100 

(Attn:  Mr.  Charles  Harvey)  .......... ........ 

Bendix  Products  Aerospace  Division,  South  Bend,  IN  46619 
(Attn:  Mr.  R.  V*  Cervelli)  ...a. 

Boeing  Aerospace  Company,  P.0.  Box  3999,  Seattle,  KA  98124 
(Attn:  Code  206,  ft?.  R.  1.  Horton) 

Boeing  Company,  Renton,  Washington  98055 

(Attn:  Dr.  R.  June) 

Boeing  Caspany,  Vertol  Divlcien,  Phils. , FA  19142 
(Attn:  Mr.  R,  L*  Pinckney,  Hr.  D.  Hoffstedt)  • 

Boeing  Company,  Wichita,  CS  67210 
(Attn:  Mr.  V.  Fcoeau/HS  16-39) 

Cabot  Corporation,  Billerica  Research  Canter,  Billerica,  HA 
01821 


Drew!  University,  folia. , PA  19104 

(Attn:  Dr.  P.  C.  Chou)  

E.L  DuPont  Company,  Wilmington,  DC  19896 

(Attn:  Dr.  Carl  Zveben)  Bldy.  262/Room  316  ..•••••••• 

Fairchild  Industries,  Hager? town,  KD  21740  * . * 

(Attn:  Hr.  D.  Ruck).  

Georgia  Institute  of  Technology,  Atlanta,  GA 
(Attn:  Prof.  W.  H.  Horton) 

General  Dynamics /Cotrvalr,  San  Diego,  CA  92138 

(Attn:  Hr.  D.  8.  Dunbar,  V.  C.  Scbeck)  ••.•••••»*•• 
General  Dycsaics,  Fort  Worth,  TX  76101 

(Attn:  Mr.  P.  D.  Stuckey,  Dept.  23,  Kail  Zone  P-46).  . . , . , 
General  Electric  Ccspany,  Phils.,  PA  19161 

(Attn:  Mr*  La  McCreight)  . , 

Great  Lakes  Carbon  Corp.,  H.Y,  , NT  1001  V 

(Attn:  Mr.  W.  R.  Bans,  Karkey  Development)  ..••••• 

Grumman  Aerospace  Corporation,  Bathpage,  L.X.,  HI  11714 
(Attn:  Mr.  R.  Hadcock,  Mr.  S,  Das  tin). 
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Non-Goveroaent  Agencies  (Cont.) 


Hercules  Powder  Company,  Inc,,  Cumberland,  HD  21501 

(Attn:  Mr.  D.  Eug) 1 

H.  I.  Thompson  Fiber  Gloss  Company,  Gardena,  CA  90249 

(Aten:  Mr.  N.  Hyets)  ........  1 

ITT  Research  Institute,  Chicago,  IL  60616 

(Attn:  Mr.  K.  Hofar)  • *••••••••••••••  1 

J.  P.  Stevens  & Co.,  Inc.,  N.Y. , MY  10036 

(Attn;  Mr.  H.  I.  Shu lock).  1 

Caman  Aircraft  Corporation,  Bloomfield,  CT  06002 

(Attn:  • Tech.  library).  ...r.  ......  ••••«.•••  1 

Lehigh  University,  Bethlehem,  PA  18015 

(Attn:  Dr.  G.  C.  Sih) ••••.•••••••«•  1 

Lockheed-Callfornia  Company,  Burbank,  CA  91520 

(Attn:  Mr.  E„  K.  Walker,  R.  L.  Vaughn) 2 

Lockheed -Georgia  Company,  Marietta,  GA 

(Attn:  Advanced  Composites  Information  Center,  Dept.  72-14, 

Zone  42).  ••••••••»••  1 

LTV  Aerospace  Corporation,  Dallas,  TX  75222 

(Attn:  Mr.  0.  E.  Dhonau/2-53442 , C.  R.  Foreman). • 2 

Martin  Company,  Baltimore,  MD  21203 

(Attn:  Mr.  J.  E.  Pawken) «••»•••  1 

Materials  Sciences  Corp. , Blue  Bell,  PA  19422  .•••••«•••  1 

McDonnell  Douglas  Corporation,  St.  Louis , M0  63166 

(Attn:  Mr.  R.  C.  Goran,  O.  B.  McBee,  C.  Stenberg).  •••»••  3 

McDonnell  Douglas  Corporation,  Long  Beach,  CA  90801 

(Attn:  H.  C.  Schjulderup,  G.  Lehman)  .«•••••••••••  2 

Minnesota  Mining  and  Manufacturing  Company,  St.  Paul,  MM  55104 

(Attn:  Mr.  W.  Davis)  1 

Northrop  Aircraft  Corp.,  Horair  Dlv. , Hawthorne,  CA  90250 
< (Attn:  Mr.  R.  D.  Hayes,  J.  V.  Noyes,  R.  C.  Isemann).  • • « • . 3 

Rockwell  International,  Columbus,  OH  43216 

(Attn:  Mr.  0.  G.  Acker,  K.  Clayton). ••••••  2 

Rockwell  International,  Los  Angeles,  CA  90053 

(Attn:  Dr.  L.  Packman)  • • • 1 

Rockwell  International,  TVlsa,  OK  74151 

(Attn:  Mr.  E.  Sanders,  Mr.  J.  H.  Powell)  .*...»•••••  2 
Owens  Corning  Fiberglass,  Granville,  OH  43023 

(Attn:  Mr.  D.  Mattes).  • 1 

Rohr  Corporation,  Riverside,  CA  92503 

(Attn;  Dr.  F.  Riel  atid  Mr.  R.  Elkin)  • •••.«•••  2 

Ryan  Aeronautical  Company,  San  Diego,  CA  92112 

(Attn:  Mr.  R.  Long).  1 

Sikorsky  Aircraft  * Stratford,  CT  06497 

(Attn:  Kr.  J.  Ray) ••••* . 1 

University  of  Oklahoma,  Norman,  OK  93069 


Union  Carbide  Corporation,  Cleveland,  OH  44101 

(Attn:  Dr.  B.  F.  Volk)  1 
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